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An  issue  of  paramount  importance  in  the  study  of  large  scale 
systems  is  that  of  model  simplification  or  reduced  order  modeling.  The 
sheer  size  on  the  one  hand  and  the  richness  and  complexity  of  phenomena 
on  the  other  make  the  use  .of  detailed  models  in  the  analysis  and  control 
of  large  systems  impractical  if  not  impossible.  A  good  analyst  or  designer 
knows  that  a  model  should  encompass  only  the  "relevant"  behavior  of  the 
system  and  should  not  be  cluttered  with  unnecessary  detail.  Although  this 
may  sometimes  be  accomplished  by  employing  parsimonious  models  for  the 
components  of  the  system,  there  are  cases  where  further  simplification  is 
needed  to  make  the  model  manageable  both  computationally  and  conceptually. 

A  characteristic  example  arises  in  stability  studies  of  interconnected 
power  systems  where  the  use  of  the  crudest  model  for  each  generator  (the 
so-called  electromechanical  model)  results  in  hundreds  or  even  thousands 
of  state  variables.  It  is  thus  desirable  to  have  systematic  model  order 
reduction  methods  for  which  the  approximation  involved  can  be  estimated. 

Singular  perturbations  is  a  well  documented  [1-4]  method  for 
reduced  order  analysis  and  design,  in  which  dynamic  phenomena  of  widely 
different  speeds  are  treated  separately.  In  the  short  run  the  slow 
dynamics  are  essentially  constant  and  the  focus  is  on  the  fast  ones. 

In  the  long  run  the  fast  dynamics  settle  to  their  "quasi-steady-state" 
and  the  focus  is  on  the  slow  dynamics.  This  time- scale  thinking  is  common 


in  diverse  engineering  fields  [5-7].  If  a  small  parameter  e  representing 
the  speed  ratio  of  slow  and  fast  dynamics  can  be  identified  this  intu¬ 
itively  appealing  idea  leads  to  asymptotic  analysis.  Most  of  the 
literature  [8-10]  is  devoted  to  systems  of  the  form 


y  ■  f(y,2,e) 


y(0)  =  y„ 


ez  =  g(y,z, e) 


z(0)  =  z 


where  e  multiplies  the  z-derivatives,  y  is  a  v-vector  and  z  is  a  p-vector. 
Formally  setting  e=0  in  (1.1),  solving 


0  =  g(y,z,0) 


for  z,  z  =  Y(y)  and  substituting  into  (1.1) 


=  f(y,  Y(y)  >0)-f (y) ,  y(0)  *  yc 


we  obtair  he  slow  reduced  model.  Writing  system  (1.1)  in  the  "stretched" 
time  variable  t  =  ~  and  setting  e=0  we  obtain  the  fast  reduced  system  (or 
associated  system  or  boundary  layer  system) 


TT  ■  8(y  .  z  +  z,  0)  z(0)  *  z-z 

UT  O  O  0  0 


where  zq“z(0).  Variables  y,z,e,t  are  restricted  to  lie  in  a  domain 
D  :  |jy-y(t)||  <  r,  ||z-z(t)||  <  r,  0  <  e  <  eQ,  0  <  t  <  T,  where  r  >  l'yo~y (0) !| 
The  following  theorem  relates  the  solutions  of  (1.1)  with  the  solutions  of 
(1.3)-(1.4). 


Theorem  1.1  [35]  Let  the  following  conditions  be  satisfied. 

HI.  f,  df/dy,  df/dz,  S»  bs/by,  bg/bz  are  of  class  C°  in  D. 

H2.  The  solution  'z(j )  of  (1.4)  exists  on  t€  [O,09  ],  is  unique,  and  is 
asymptotically  stable  with  respect  to  ”z=0. 

H3.  The  solution  y(t)  of  the  reduced  system  (1.3)  exists  and  is  unique 
on  tg[0,T] . 

H4.  The  real  parts  of  the  eigenvalues  of  the  Jacobian  matrix 

ag/az(y,z,0)  (1.5) 

are  negative  on  [0,T],  for  "z=Y(y). 

Then  for  sufficiently  small  e,  the  full  system  (1.1)  has  a  unique 
solution  y(t,£)  on  t£  [0,T]  satisfying  the  initial  condition  y(0,e)=yo> 
z(0,e)=zq.  Furthermore, 

lim  y(t,e)=y(t ,  on[0,T]  (1.6) 

e-»  0 

lim  z(t,e)  =  z(t)+z(-)  on  [0,T]  (1.7) 

e-0  6 

where  the  limits  in  (1.6),  (1.7)  are  uniform  in  t  on  [0,T]. 

From  (1.6),  (1.7)  the  response  of  y  in  (1.1)  is  approximated,  to 
0(e),  hy  the  response  of  the  slow  system  (1.3)  whereas  the  response  of  z 
is  approximated  by  a  boundary  layer  zf(~)  superimposed  on  the  quasi-steady- 
state  z(t)=Y(y(t)). 

An  extensive  literature  dealing  with  system  (1.1)  and  the 
corresponding  controlled  system  includes  results  on  stability  [20,65,66], 
linear  [4,67]  and  nonlinear  [3]  regulator  design,  controllability  properties 
[37]  and  time-optimal  control  [68] , filtering  and  smoothing  [69].  A  basic 


assumption  in  these  references  is  that  the  Jacobian  matrix  3g/dz(y,z,0) 
is  nonsingular.  When  this  happens  we  say  that  time  scales  in  (1.1)  are 
explicit ,  that  is,  they  coincide  with  the  decomposition  of  the  state  vector 
into  y  and  z.  When,  however,  Sg/^z  is  singular  time  scales  in  (1.1)  are 
nonexplicit,  that  is,  all  states  may  be  mixed  having  fast  and  slow  parts. 

Some  authors  treat  such  cases  as  "singular-singularly  perturbed"  systems 
[11-13]  or  "generalized  singularly  perturbed"  systems  [14,15]. 

In  this  thesis  we  take  an  alternate  route.  We  recognize  that  in 
a  wide  class  of  systems  singularity  of  gg/^z  Is  due  to  the  selection  of 
state  variables;  hence  a  nonsingular  transformation  removes  the  singularity 
of  dg/dz,  and  defines  new  states  in  which  time-scales  are  explicit.  This 
approach  has  two  advantages.  First  it  puts  the  system  into  a  form  in  which 
the  results  alluded  to  before  can  be  applied.  Second,  from  the  transformed 
system  we  can  easily  define  fast  and  slow  reduced  systems  describing  the 
system  behavior  in  the  short  run  and  in  the  long  run.  Following  this 
approach  we  establish  a  relation  between  weak  connections  and  time  scales 
in  a  class  of  interconnected  systems.  Separation  of  time  scales  in  such 
systems  leads  to  a  physical  decomposition  into  a  slow  core  and  a  number  of 
weakly  coupled  fast  subsystems.  The  results  are  further  specialized  to 
structured  interconnected  systems  such  as  power  systems  and  other  dynamic 
networks. 

1.2  Chapter  Preview 

Chapter  2  starts  with  a  simple  RC  example  pointing  the  relationship 


between  equilibrium  and  conservation  properties  on  the  one  hand  and  time 
scales  on  the  other  [16].  These  properties  are  used  in  the  construction  of 


a  transformation  that  makes  the  time  scales  in  Linear  Time  Invariant  (LTI) 
systems  explicit.  Next  a  multi-time-scale  system  is  viewed  as  a  succession 
of  two -time -scale  ones.  Starting  from  the  fastest  time  scale  we  proceed 
to  the  slower  ones,  using  at  each  step  the  transformation  that  makes  time- 
scales  explicit.  This  procedure  defines  a  sequence  of  "nested"  reduced 
order  models.  The  transformation  separating  the  time  scales  is  then 
generalized  to  LTI  systems  with  inputs. 

In  Chapter  3  the  equilibrium  and  conservation  reasoning  is  extended 
to  nonlinear  systems  leading  to  a  transformation  that  makes  time  scales 
explicit  in  models  of  the  form 

ex  =  h(x,e) .  (1.8) 

It  is  then  shown  that  the  results  of  [17-19]  on  high  gain  feedback  and 
disturbance  decoupling  generalize  to  a  class  of  nonlinear  systems  for 
which  the  controls  enter  linearly  but  the  output  map  and  feedback  law  are 
nonlinear.  We  next  turn  to  interconnected  systems  made  of  systems  with 
equilibrium  manifolds  and  show  that  weak  connections  give  rise  to  two-time¬ 
scale  behavior.  A  decomposition  of  interconnected  systems  into  a  slow  core 
and  fast  local  systems  leads  to  decentralized  stability  criteria  based  on 
the  results  of  [20] . 

Chapter  4  deals  with  time  scales,  coherency  and  aggregation  in 
nonlinear  dynamic  networks.  Coherency  based  aggregation  [21-23],  a  common 
procedure  for  order  reduction  in  power  systems,  is  given  theoretical 
foundations  for  nonlinear  electromechanical  models,  thus  extending  the 
results  of  [24-29] .  It  is  shown  that  linear  physical  laws  result  in 
linear  time-scale  separating  transformation  even  when  some  components  of 


6 

the  network  are  nonlinear.  A  five-machine  power  system  example  illustrates 
the  proposed  reduced-order  modeling  and  verifies  its  validity. 

Extensions  in  several  directions  and  possible  uses  of  the 
decomposition  in  direct  stability  analysis  are  discussed  in  Chapter  5. 


CHAPTER  2 


MODELING  OF  TWO -TIME-SCALE  SYSTEMS 

2 . 1  Introduction 

When  the  model  of  a  real  system  with  the  two-time-scale  property 
is  expressed  in  terms  of  physical  variables  it  often  fails  to  be  in  the 
form  (1.1).  An  important  requirement  in  (1.1)  is  that  dg/^z  be  nonsingular 
along  li(t).  When  this  condition  is  violated  the  model  is  said  to  be  non¬ 
explicit  and  the  conclusions  in  Chapter  1  have  to  be  modified. 

Some  authors  [11-13]  treat  nonexplicit  models  as  "singular- 
singularly  perturbed"  systems.  Instead  we  approach  them  from  the  modeling 
point  of  view  recognizing  that  the  singularity  of  dg/dz  is  due  to  the  choice 
of  state  variables.  We  show  that  equilibrium  [16]  and  conservation  properties 
provide  a  coordinate-free  characterization  of  singular  perturbations.  These 
properties  are  used  in  the  construction  of  a  transformation  leading  to  the 
explicit  model  (1.1)  with  the  slow  part  of  z(t)  being  0(e).  The  discussion 
in  this  chapter  is  restricted  to  Linear  Time  Invariant  (LTI)  systems. 

Extension  of  the  basic  ideas  to  nonlinear  systems  and  applications  to 
interconnected  systems,  high  gain  feedback  and  dynamic  networks  appear  in 
Chapters  3  and  4. 

In  Section  2.2,  a  simple  physical  system  is  used  to  motivate  the 
discussion  and  indicate  the  relation  between  time  scales  on  the  one  hand 
and  equilibrium  and  conservation  properties  on  the  other.  In  Section  2.3,  the 
relation  is  established  for  LTI  systems  and  a  transformation  is  constructed 
that  transforms  a  nonexplicit  singularly  perturbed  model  to  the  explicit 


model  (1.1).  In  Section  2.4,  multi-time-scale  systems  are  treated  as  a 
succession  of  two-time-scale  systems  and  a  sequence  of  nested  reduced  order 
models  is  defined.  Section  2.5  generalizes  the  results  of  Section  2.3  to 


systems  with  inputs  and  Section  2.6  deals  with  some  structured  nonexplicit 
models. 

2.2  Equilibrium  and  Conservation  Properties 

Although  nonexplicit  singular  perturbations  occur  in  as  simple 
systems  as  RC-circuits  they  have  not  attracted  much  attention.  In  contrast, 
explicit  perturbations  have  been  investigated  for  networks  with  "parasitic" 
inductances  and  capacitances  [30,31].  When  such  parasitics  are  expressed 
as  multiples  of  e  capacitor  voltages  and  inductor  currents  are  used  as 
state  variables,  the  circuit  model  is  in  the  explicit  form  (1.1).  A  simple 
illustration  is  the  RC-circuit  of  Fig.  2.1a  with  state  equations 

(RjCp  dx^/dt^  -  -xL  +  x2 

(R^)  dx2/dtd  -  Xj_  -  [1  +  (Rj/Rj)]  x2  +  <W  x3  (2-1) 

(R1C3)  dx3/dtd  -  (Rl/R2)  x2  -  [(Rj/R2)  +  (R1/R3)]  x3 

where  the  capacitor  voltages  were  chosen  as  states  and  td  is  dimensional 
time.  Suppose  that  C2  and  C3  are  "parasitic,"  say  C2=C3=eC^  and  that 
all  the  resistors  are  of  the  same  order  of  magnitude.  Recognizing  RjC^ 
as  a  typical  large  time  constant  and  defining  the  slow  dimensionless  time 
t*t./(R,C,)  [32,33],  (2.1)  becomes 


s 


ed^/dt  -  xL  -  [1  +  (Rx/R2)]  x2  +  (R^)  x3  (2. 

edx3/dt  -  (Rl/R2)  x2  -  [(R^)  +.  (Rj/R^I  x3 

with  e  multiplying  the  derivatives  of  X2  and  Xy  Thus  x^  appears  as  the 
y-variable  and  x2>x3  as  t*ie  z-variables  of  (1.1)  and  the  model  is 
explicit  because  the  two-by-two  matrix  of  x2»x3  *-s  nonsingular.  The  slow 
reduced  model  (1.2)  represents  the  circuit  with  parasitic  capacitors  C2 
and  C3  opened  (Fig.  2.1b),  whereas  in  the  fast  reduced  model  (1.3)  the 
large  capacitor  C3  is  shortened  [30]  (Fig.  2.1c). 

In  the  same  circuit  nonexplicit  singular  perturbations  occur 
when  all  capacitors  are  of  the  same  order  of  magnitude,  say  C^=C2“C3=C, 
but  the  resistors  are  not.  For  example,  if  R^  and  R2  are  small  and  R3 


is  large,  say 


R.  ■=  r,  R  =  r/2,  R,  =  R,  r/R  =  e 


typical  large  and  small  time  constants  are  RC  and  rC,  respectively,  and 
in  the  dimensionless  time  variables 


c  =  rc’  T  a  7c’  7  =  € 


the  circuit  is  described  by 


e  (dx/dt)  ■  dx/dT 


"-1  1  0  “ 

* 

1-3  2  x 


A(e)  x. 


Note  that  e  multiplies  all  the  derivatives  in  the  slow  time-scale  t  and 

thus  there  are  no  explicit  slow  y- variables  in  the  system.  If  that 

is  A(0),  were  nonsingular,  no  slow  phenomenon  would  exist  in  (2.5)  and 

the  system  would  not  possess  the  two-time-scale  property.  However,  A(0) 

is  singular  indicating  the  existence  of  a  ''hidden"  slow  phenomenon.  To 

T 

see  this  assume  that  x(0)  =  [1  1  1]  in  (2.5).  Then  the  slow- time 

derivatives  dx/dt  remain  finite  when  e-*  0  suggesting  that  (2.5)  is  a 
two-time-scale  system.  Physically  the  slow  phenomenon  is  the  discharge 
of  the  capacitors  through  the  large  "leakage"  resistor  R^.  Neglecting 
this  "leakage,"  Fig.  2.1b,  makes  the  slow  phenomenon  infinitely  slow, 
that  is, constant  and  corresponds  to  setting  e=0  in  the  T-model  of  (2.5) 


2  x  =  A(0)x. 


Since  A(0)  is  singular  the  equation 


A(0)  x  =  0 


has  an  infinite  number  of  solutions  given  by 


x  =  a  [1  1  1] 


where  a  is  any  real  number,  that  is, (2.6)  has  a  continuum  of  equilibrium 
points.  This  can  be  seen  from  the  circuit  of  Fig.  2.1b  where  any  x  such 


(2.9 


x!  "  x2  “  °»  x3  ‘  x2  “  0 


is  an  equilibrium  point.  The  line  represented  by  (2. 8)- (2. 9)  will  be 
denoted  by  S. 

Kirchhoff’s  current  law  (KCL)  applied  to  the  ground  node  of 
Fig.  2.1b  gives  the  dual  property  that  is  the  conservation  of  total  charge 
for  all  t. 


Cj.*!.  (t>  +  C2X2(t)  +  C3x3(t)  =  Clxl(0)  +  C2X2<°>  +  C3X3(0)  (2*10 

which  means  that  every  trajectory  x(t)  of  (2.6)  is  confined  to  a  plane  F 
passing  through  the  initial  point  x(0)  orthogonal  to  the  vector 
[C^  C2  C3]  .  The  quantity  in  (2.10),  constant  when  e=0,  becomes  slowly 
varying  when  e  >  0,  that  is  when  the  "leakage"  is  introduced.  A 
circuit  describing  this  slow  phenomenon  is  given  in  Fig.  2.1c  and  will 
be  derived  in  the  next  section. 

From  the  above  discussion  we  conclude  that  the  trajectories 
x(t)  of  the  original  system  2.5  consist  of  two  distinct  parts.  First 
in  a  "boundary  layer"  near  plane  F  the  state  x(t)  rapidly  approaches 
line  S.  Then,  from  a  neighborhood  of  the  intersection  of  plane  F  with 
line  S,  x(t)  continues  to  slowly  "slide"  along  line  S.  The  geometry  of 
this  situation  is  sketched  in  Fig.  2.2.  Note  that  the  behavior  of  x(t) 
is  similar  to  that  of  the  explicit  model,  that  is,  a  fast  transient  is 
followed  by  a  slow  motion  close  to  a  line  of  "quasi-equilibria"  S.  The 
basic  difference  is  that  in  the  explicit  model  (2.2)  the  plane  F  close 
to  which  the  boundary  layer  occurs,  is  orthogonal  to  axis  x.  (Fig.  2.3). 
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Examination  of  Fig.  (2.2) -(2. 3)  indicates  that  in  nonexplicit  models  fast 
dynamics  are  observed  by  all  states  whereas  in  explicit  models  they  are 
only  weakly  observed  by  some  states  (the  y-variables) . 

This  example  indicates  that  the  time  scales  of  the  original 
system  (2.5)  are  related  to  the  equilibrium  and  conservation  properties 
of  the  auxiliary  system  (2.6)  in  T-scale.  These  properties  are  coordinate 
free  and  characterize  all  two-time-scale  systems  reducible  to  the  explicit 
model.  In  the  next  section  they  will  serve  for  a  choice  of  coordinates  in 
which  the  time  scales  are  explicit. 

2.3  Nonexplicit  Singularly  Perturbed  Systems 

The  discussion  of  the  previous  section  will  now  be  generalized 
to  the  system 


edx/dt  =  dx/dr  =  A(e)  x  (2.11) 

where  x6R.n,  A(e)  is  a  time  invariant  nxn  matrix  depending  on  e,  and  t, 
r  the  slow  and  fast  time  variables,  respectively.  The  following  is 
assumed  about  A(e). 

Assumption  2.1  A(e)  can  be  written  as 


A(e)  -  Aq  +  e  A^(e)  (2.12) 

with  A^(e)  bounded  at  e=0  and  A(0)=Ao  satisfying 

R(A0)  ©  7KAq)  -  Rn  W  (2.13) 


(*) 

In  the  terminology  of  [13],  (2. 13) - (2 . 14)  is  equivalent  to  ind 
A  *1.  In  [13]  it  is  shown  that  this  condition  is  necessary  for  lim  x(t)  to 
exist.  e-  0 
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where  R(A  )  is  the  range  space  of  A  ,  ??(A  )  is  the  null  space  of  A  ,  and 
o  o  o  o 

@  denotes  the  direct  sum  of  two  spaces  [34].  The  dimensions  of  ft(Ao), 
7KAo)  are 


dim  R(Aq)  =  p  >  1,  dim  7[(Aq)  =  v  >  1,  p  +  v  =  n. 


(2.14) 


Equation  (2.13)  is  equivalent  to  saying  that  Aq  has  a  complete 
set  of  eigenvectors  corresponding  to  its  zero  eigenvalues,  which  in  turn 
is  equivalent  to  the  following:  ft(AQ)  is  the  invariant  space  (eigenspace) 
of  Aq  corresponding  to  the  nonzero  eigenvalues,  and  7}( Aq)  is  the  invariant 
space  (eigenspace)  corresponding  to  the  zero  eigenvalues. 

To  study  the  time-scale  behavior  of  (2.11)  the  auxiliary  system 


dx/dt  *  A  x 
o 


(2.15) 


is  defined  with  A  as  in  (2.12).  By  assumption,  (2.15)  has  a  v-dimensional 
o 

/*\ 

equilibrium  manifold v  '  (i.e.  ^(Aq))  S  consisting  of  all  x  such  that 


A  x  =  0. 
o 


(2.16) 


If  W  is  a  pxn  matrix,  rank  W=p,  whose  rows  span  the  row  space  of  A,  then 


Wx  =  0 


V  x^S 


(2.17) 


To  see  the  conservation  property  of  (2.15)  we  note  that  if  V 
is  a  \>xn  matrix,  rank  V  =  v  whose  rows  span  the  left  null  space  of  Aq, 
i.e.  VAq  =  0,  then 


'  'Although  S  is  presently,  simply  a  subspace,  we  call  it  manifold 
in  anticipation  of  the  nonlinear  extension  in  Chapter  3. 
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V (dx/dr)  =  V  Aq  x  =  0  V  t,  V  x(0)€Rn.  (2.18) 

Thus,  the  v-dimensional  quantity  Vx  is  constant  along  the  trajectories 
of  (2.15), 

V  x(t)  =  V  x(0),  Vx(0)€Rn.  (2.19) 

This  means  that  for  each  value  of  V  x(0)  the  trajectory  of  (2.15)  is 
confined  to  a  linear  manifold ^  defined  by  (2.19). 

This  linear  manifold  called  dynamic  manifold  F,  is  orthogonal 
to  the  rows  of  V  and  contains  the  initial  point  x(0).  The  orthogonality 
between  the  left  null  space  and  the  range  space  of  a  matrix  [34]  implies 
that  F  is  a  translate  of  R(Aq). 

The  above  discussion  has  established  equilibrium  (Eq.  (2.17)) 
and  conservation  (Eq.  (2.19))  properties  analogous  to  the  ones  of  the 
RC-circuit  of  the  previous  section  (Eq.  (2.9)  and  (2.10)).  The  behavior 
of  the  trajectories  is  still  the  one  depicted  in  Fig.  2.2  with  S  and  F 
defined  by  (2.17)  and  (2.19),  respectively.  We  are  now  ready  to  define 
a  new  set  of  coordinates  in  which  the  time  scales  are  explicit. 

Theorem  2.2  Under  Assumption  (2.1)  the  change  of  coordinates 

y  =  Vx,  z  -  Wx  (2.20) 

transforms  (2.11)  into  the  explicit  model  (1.1)  with  'z(t)  =  0. 


v  ''A  linear  manifold  of  dimension  r  is  a  translation  of  an 
r-dimensional  subspace. 
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Proof:  Since  the  rows  of  V,W  form  bases  for  the  left  null  and  row  spaces 


of  Aq,  respectively,  the  transformation 


(2.21) 


defined  by  (2.20)  has  inverse 


T  -  [P  Q] 


(2.22) 


where  the  columns  of  P,Q  form  bases  of  72(Aq),  R(Aq).  Hence, 


T  (rf  +  a1(€))t 


WAP  W  A  Q 
o  o 


-1  1  V  Ao  P  V  Ao  Q  +  V  Ai<e>P  v  al(OQ 


WAx(e)P  WA1(e)Q 


VA^e)?  V  Al(e)Q 


W  A  Q 

WAl(e)P  ~  +  WA1(£)Q 


(2.23) 


dy/dt  =  A1]L(e)  y  +  AL2(e)  z 


edz/dt  =  e  A2L(e)  y  +  A22(e)  z 


(2.24) 


where  An(e)  A  v  A1(e)P,  A12(s)  A  v  A1(e)Q,  A21(e)  ^  W  A^e)?  and 


A,,(e)  =  W  A  Q  +  e  W  A.(e)Q. 


(2.25) 


To  show  that  (2.24)  is  explicit  model  we  need  to  show  that  A22(0)  is 
nonsingular.  Notice  that  Assumption  (2.1)  implies  that  R(Aq)  is  the 
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eigenspace  of  the  nonzero  eigenvalues  of  AQ.  Hence  there  is  a  pxp 

nonsingular  matrix  G  whose  eigenvalues  are  the  nonzero  eigenvalues 

of  A  such  that 
o 


A  Q  =  QG. 
o 


(2.26) 


The  last  relation  implies 


A„(0)  -  WAQ-VQG-G 

22  O 


(2.27) 


which  is  nonsingular. 
Remark:  Writing 


=  T 


1  r y_ 
2  - 


=  Py  +  Qz 


(2.28) 


we  see  that  y,z  are  the  representations  of  x  with  respect  to  bases  P,Q 
of  72(Aq),  R(Aq),  respectively.  Another  way  to  view  (2.28)  is  that  Py 
is  the  projection  of  x  on  7?(Aq)  along  R(Aq)  and  hence  y  is  the  representa¬ 
tion  of  this  projection  with  respect  to  basis  P.  A  similar  interpretation 
holds  for  Qz  and  z. 

We  now  illustrate  the  application  of  Theorem  2.2  with  the  RC-circuit 
of  Fig.  2.1  in  which  the  time  scales  are  due  to  large  and  small  resistors  as 
in  (2.3).  The  auxiliary  system  is  given  in  (2.6)  and  the  equilibrium  and 
dynamic  manifolds  are  defined  by  (2.9)  and  (2.10),  respectively.  A  choice 
of  coordinates  according  to  these  equations  is 


y  ■  +  ^2*2  ^3x3^^a 


(2.29) 


where  the  division  by 


Cl  + 


C„  +  C, 


(2 . 31 


in  (2.29)  retains  the  physical,  meaning  of  y  as  a  voltage  variable.  In  the 
new  coordinates  the  circuit  is  described  by 

dy/dt  =  -  (C/Ca)  y  +  (C2/Ca)  Zj.  -  (2C2/C2)  ^ 

e(dZ]L/dt)  =  -  2zl  -  2z2  (2.32; 

e(dz2/dt)  =  -  ey  -  (1  -  e(C/Ca))  -  (4  +  e(2C/Ca))  z2- 

As  stated  in  Theorem  2.2  the  z-equations  in  (2.32)  give  "z(t)  =  0.  Hence 
the  slow  reduced  model  is 


dy/dt  =  -  (C/C  )y 

a 


(2.33; 


represented  by  the  circuit  in  Fig.  2.1c  and  the  fast  reduced  model  is 


dzL/dT  -  -  2zL  -  2z2 
dz2/dx  =  -  z'1  -  4z2 


(2.34] 


represented  by  the  circuit  in  Fig.  2.1b  where  the  voltages  with  respect  to 
the  "reference"  node  2  are  used  as  states. 

It  is  interesting  to  note  the  physical  interpretation  of  the 
new  variables.  The  slow  variable  y  is  proportional  to  the  sum  of  the 
charges  on  the  three  capacitors  and  can  be  considered  the  voltage  on  the 
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“aggregate"  capacitor  C  (Eq.  (2.31),  Fig.  2.1c).  Application  of  Kirchhoff’s 
current  law  to  the  ground  node  of  Fig.  2.1a  shows  that  the  time  derivative  of 
y  is  proportional  to  the  current  in  R^>  Since  is  large  dy/dt  is  small 
and  y  qualifies  as  a  slow  variable.  The  fast  variable  z^  equals  the  voltage 
across  R^  which  due  to  the  smallness  of  R^  diminishes  quickly  to  values 
close  to  zero.  Hence  z^  qualifies  as  a  fast  variable.  A  similar  interpreta¬ 
tion  holds  for  z2‘  In  Chapter  4  we  show  that  this  selection  of  slow  and 
fast  variables  is  good  for  a  wide  class  of  nonlinear  dynamic  networks. 

2.4  Nested  Reduced  Order  Models 

In  the  previous  section  we  showed  how  equilibrium  and  conservation 
properties  are  used  to  transform  a  nonexplicit  singularly  perturbed  model  to 
an  explicit  one.  Writing  the  explicit  model  (2.24)  in  the  fast  time  t  and 
letting  e-*  0  we  obtain  the  fast  reduced  model 

g=A22(0)z  z(0)  =  z(0).  (2.35) 

Similarly  writing  the  model  in  the  slow  time  t  and  letting  e-*  0  we  obtain 
the  slow  reduced  model 

&  =  Au(0)7  y(0)  =  y(0) .  (2.36) 

If  the  eigenvalues  of  A22 (°)  have  negative  real  parts  the  responses  of 
(2.35)-(2.36)  are  0(e)  approximations  of  the  response  of  (2.24)  over  bounded 
intervals  [8-10].  If,  in  addition,  the  eigenvalues  of  the  slow  system 
matrix  A^(0)  have  negative  real  parts  the  approximation  is  valid  over 
unbounded  intervals  [36]. 


l: 
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It  may  happen,  however,  that  some  eigenvalues  of  A^(0)  are  zero. 

In  these  cases  the  approximation  is  not  valid  over  unbounded  intervals 
since  the  response  of  (2.36)  tends  to  a  nonzero  constant  whereas  the  response 
of  (2.24)  tends  to  zero.  Treatment  of  (2.24)  as  a  two-time-scale  system 
is  inadequate  because  the  system  may  have  more  than  two  time  scales.  Simple 
expansions  of  the  eigenvalues  and  eigenvectors,  for  example,  show  that  if 
Aj^  is  singular  the  matrix 


eAll 

£A12 

eA21 

<N 

CM 

< 

which  is  the  matrix  of  (2.24)  in  the  fast  time  scale  with  A^,  A A21* 

2 

independent  of  e,  has  0(e  )  eigenvalues  in  addition  to  0(1)  and  0(e) 

A  4 

ones. 

Instead  of  treating  the  original  system  (2.11)  as  a  multi-time- 

scale  one,  we  prefer  to  deal  with  only  two  time  scales  at  a  time.  That  is, 

starting  from  the  fastest  time  t^-T  we  consider  the  system  operating  in 

scale  t^  and  t2=eti  only.  Viewed  from  t^  speeds  are  0(1)  and  the  rest  are 

2 

o (1) ;  we  do  not  specify  whether  they  are  0(e),  0(e  ),  etc.  Changing  time 
scales  to  the  slower  t^,  some  speeds  are  o(~)  and,  in  time-scale  t^,  are 
assumed  to  reach  their  quasi-steady-state  instantaneously.  The  rest  of  the 
system  is  again  treated  as  two-time-scale  with  t^  the  fast  time. 

To  make  this  idea  precise  we  employ  the  block  diagonalizing 


transformation  [37,38] 
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(2.37) 


where  L,H  satisfy 

V‘>  +  £  L  AH<£>  -  A22<£>  1  L  A12<s>  L  '  0  <2'38) 

and 

A12(e)  -  H  (A22(e)  +  e2  L  Al2(*))  +  e(An(e)  -  e  Al2(e)  L)  =  0.  (2.39) 

In  the  §,T|  coordinated  (2.24)  becomes 

|  *  (Au(e)  •  e  A12(E)  L)5  (2.40) 

e  -  (A22(e)  +  e2  L  A12(«))T|  .  (2.41) 

Since  A22(0)  is  nonsingular  (2.35)  is  a  regular  perturbation  of  (2.41)  and 
can  be  used  as  a  fast  reduced  order  model.  However,  if  A^(0)  is  singular 
satisfying  Assumption  2.1,  (2.40)  is  nonexplicit  singularly  perturbed  model 
in  the  form  (2.11).  Hence,  arguing  as  in  Section  2.3,  we  apply  transforma¬ 
tion  (2.20)  to  define  slow  and  fast  variables  in  time  scale  t2>  The  process 
can  be  repeated  until  all  time  scales  are  "peeled  off"  defining  nested 


reduced  order  models. 


Although  some  control  design  work  has  been  done  for  generalized 
singularly  perturbed  systems  [13,14],  most  of  the  literature  [2-4,39-40] 


deals  with  the  explicit  model 

y  -  Au  y  +  A12  z  +  Bj  u 
Ez  -  A2l  y  +  A22  z  +  B2  u 


(2.42) 


where  A22  is  nonsingular.  Note  that  in  the  y-equations  the  gain  of  the 
control  u  is  0(1)  whereas  in  the  z-equations  it  is  0(^-).  We  are  interested 
in  conditions  under  which  the  nonexplicit  model 

edi "  fr =  [Ao +  eVe)lx  +  ^Bo +  u  (2*43> 

where  Aq  satisfies  Assumption  2.1  and  B(e)  is  differentiable  at  e=0,  can 
be  transformed  to  the  explicit  model  (2.42). 

Assumption  2.3  Let  V  be  a  vxn  matrix  that  spans  the  left  null  space  of 

A  .  Then 
o 

V  B  *  0  (2.44) 

o 

that  is,  V  is  in  the  left  null  space  of  Bq. 

Corollary  2.4  Under  Assumptions  2.1,  2.3  the  transformation 

y  ■  Vx,  ,  z  *  Wx  (2.45) 

transforms  (2.43)  into  the  explicit  model  (2.42)  with  "z  =  -  A22(0)B2(Q)u 
where  B«(e)  *  WB  +  cW5(e)  and  u  is  the  slow  control. 
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Proof:  Follows  directly  from  Theorem  2.2  and  (2.44). 

Condition  (2.44)  essentially  requires  that  the  control  driving 
the  slow  variable  y  have  0(1)  gain,  as  in  the  explicit  model  (2.42). 

If  this  condition  is  not  met  the  variable  y,  slow  in  the  free  system 
(2.11),  is  subjected  to  high  gain  control  altering  the  time-scale  behavior 
of  the  system.  We  will  have  more  to  say  about  high  gain  feedback  in  the 
next  chapter.  Condition  (2.44)  is  likely  to  be  satisfied  in  well  defined 
physical  problems  as  demonstrated  by  the  following  example. 

Consider  the  transformer  of  Fig.  2.4  where  L^,  are  the  self¬ 
inductances  of  the  coils,  M  is  the  mutual  inductance  and  x^,  x^  are  the 
currents  through  the  coils.  Using  x^,  x2  as  states,  the  state  description 
of  the  system  is 


dxL/dt  =  -  (R1L2/d)x1  +  (MR2/d)x2  -  (I^/d)  v± 
dx2/dt  =  (MR1/d)x1  -  (R2L1/d)x2  +  (M/d) 

2 

where  d  =  ~  ^  ' 


In  the  case  of  an  ideal  transformer,  d  = 


M  =  0. 


(2.46) 


For 


nonideal  transformer  with  small  leakage 


lil2  ;  M 

L1L2 


=  e 


(2.47) 


where  e  is  a  small  positive  parameter.  Using  (2.47)  the  system  matrix 
of  (2.46)  becomes 


or,  substituting  */l-e  =  1  -  j  +  0(e  ) 

A(e)  *  7  (Aq  +  eA1(e)) 


where 


A 

o 


R2 

^1L2 


The  left  null  space  of  Aq 


V  =  [ll  JT^l2) 


and  the  row  space  of  Aq 

W  =  [rl 


define,  according  to  (2.45),  the  slow  variable 


y  =  Llxl  +  x2 


and  the  fast  variable 


In  the  new  coordinates  y,z  the  state  equations  become 


fj:  -  -  [1/(T1+T2)]y  +  [(T1-T2)/2(T1+T2)]z  -  ^ 


edt  £^T2  ’  TL^2(T1+T2^y 


[-1-  +  -i - § — 

lT  T  T  +T 

il  2  12 


]z 


,_L  .  U 
(t  T  3 
X1  l2 


(2.55 


where  *  L^/R^,  T2  »  L2^R2  are  t*ie  fc*me  constant:s  °f  tlie  primary  and 
secondary  R-L  circuits.  The  slow  model  is  obtained  by  setting  e=0  in 
the  second  equation  of  (2.55)  and  substituting  the  quasi-steady-state 


z  =  -  v. 


(2.56) 


into  the  first  equation  giving 


= 

dt 


l- 


1  3y 


Tl+T2 


[— ]  v 
lT1+T2j  i 


(2.57) 


The  fast  model  is  obtained  by  writing  the  second  equation  of  (2.55)  in 
the  fast  time-scale  and  setting  e=0 


dz 

dT 


,  1  1 

(=7-  +  m~)  z 

A1  l2 


(2.58] 


It  is  interesting  to  note  the  physica-1  interpretation  of  the  new  variables 
y,z  and  of  (2.54).  By  writing 


y  -  L.x  +  JlTT  x  -  L.x  +MX  +  0(e)  *  <$,,  +  +  0(e) 


(2.59 


we  see  that  y  is,  to  0(e),  the  total  flux  linkage  (j>^  +  ^12  *n  coi^ 
Aggregate  physical  quantities  such  as  total  flux  linkage,  total  charge, 
total  momentum  etc.,  are  often  slow  variables. 


Noticing  that  t^le  v°ltaSe  V£  across  the  secondary  winding 


of  the  transformer  and  using  (2.54,  (2.56)  we  obtain  as  e-*  0 


V2  -  -  (V^R^) 


■^2  ^2 
£~  =  -  (Vj^+Rj^Xj^)  — 


(2.60) 


where  are  the  number  of  turns  in  coils  1  and  2.  When  the  voltage 

drop  R.x  is  small  compared  to  v.,  which  is  usually  the  case,  (2.60) 

11  1  N2 
reduces  to  the  corresponding  relation  v^  =  -  —  v^  for  a  leakage-free 

transformer. 

Writing  M  =  «/I~£  ~  VLji^  (1  -  j)  we  see  that,  for  this 

example.  Assumption  2.3  is  satisfied. 


C 


2.6  Structured  Singularly  Perturbed  Forms 

Section  3  dealt  with  the  nonexplicit  singularly  perturbed  form 
(2,11),  in  which  all  the  states  are,  generally,  mixed.  However,  it  is 
well  known  that  a  more  structured  system  matrix  implies  that  state  x^ 
is  predominantly  slow.  In  this  section,  the  methodology  developed  in 
Section  3*  is  used  to  study  two  other  structured  forms,  the  fast  separated 
form  and  the  weak  connection  form. 

A  system  is  said  to  be  in  the  fast  separated  singularly  perturb*. 


form  if 


span  the  left  null  and  row  spaces  of  A  ,  respectively,  and  that 
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a  : 


n  i 


A  A"1 

12  22 


(2.66) 


span  the  null  and  range  spaces  of  Aq,  respectively.  Morover  they  satisfy 
(2.13). 

Corollary  2.5  If  A22  is  nonsingular,  then  (i)  the  change  of  coordinates 


y  =  xl  -  A12A22  x2,  z  =  x2 


(2.67) 


transforms  (2.61)  into  the  explicit  model  (1.1)  with  z(t)  =  0.  (ii)  The 
slow  reduced  model  of  (2.61)  is 


dt  (A11  "  A12A22  A21)y 


(2 . 68) 


and  the  fast  reduced  model  is 


dz  .  ~ 
dT  A22  Z 


(2.69) 


(iii)  The  state  x2  of  (2.61)  is  predominantly  fast  whereas  x^  is  mixed. 
Proof;  (i)  follows  directly  from  Theorem  2.2,  using  (2.64),  (2.65).  (ii) 
is  obtained  by  bearing  in  mind  that  "z(t)  =  0.  (iii)  follows  by  inverting 
(2.67) 


*1  -  y  +  A12A22  X, 


x2  =  z 


(2.70) 


Note  that  state  x2  (=z)  can  be  used  in  the  fast  reduced  model  (2.69) 
justifying  the  name  "fast  separated." 


We  now  turn  to  the  weak  connection  form  which  arises  naturally  in 


dynamic  networks  made  of  weakly  connected  "areas"  [24-29] .  A  system  is 

:(*) 


said  to  be  in  the  weak  connection  form  if 


dxl 

1 

o. 

X 

1— 

_ 1 

dt 

= 

dr 

= 

dx2 

dx^ 

_  dt  _ 

dr 

All+eAll 


eA 


12 


eA, 


21 


A22+eA22 


(2.71) 


where  x^  and  x2  are  n^  and  vectors,  A^,  A^,  A12,  A2l,  A22 


are  matrices  of  appropriate  dimensions  and  A^,  A22  are  singular  matrices 
with  a  complete  set  of  eigenvectors  corresponding  to  the  zero  eigenvalues, 
that  is,  satisfying 


n. 


n^ 


ft(An)  ©  7Z(An  )  =  R  L  , 


11' 


R(A22)  ©  7?(A22)  =  R 


(2.72) 


with 


dim  R(An)  ■  pL  >  1,  dim  >  1,  PL  +  ^  =  ^  (2.73) 

dim  R(A22)  =  p2  >  1,  dim  7KA22)  *  >  1,  p£  +  v2  =  n2  .  (2.74) 


Writing 


All+eAll 


eA 


12 


eA 


21 


A22+eA22 


A11  ° 


0  A 


22 


+  e 


A11  A12 


A21  A22 


(2.75) 


(*), 


For  convenience  we  deal  with  only  two  weakly  connected  "areas" 
but  the  ideas  are  directly  applicable  to  any  number  of  areas. 
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we  obtain 


A  = 
o 


(2.76) 


which,  because  of  (2.72),  satisfies  Assumption  2.1.  Let  V.,  span  the 
left  null  and  row  spaces  of  A^  and  Qt  span  the  null  and  range  spaces, 
respectively,  i=l,2,  satisfying 


Vi  =  Vi  =  ^i* 


i=l,2  . 


(2.77) 


Corollary  2.6  Under  assumption  (2.72),  (i)  the  change  of  coordinates 


Vi 


wlxl 


(2.78) 


transforms  (2.71)  into  the  explicit  model  (1.1)  with  z(t)  =  0.  (ii)  The 
slow  reduced  model  is 


V1  A11  P1 


V2  A21  P1 


V1  A12  P2 


V2  A22  P2 


(2.79) 


and  the  fast  reduced  model  is 


WlAllQl+eWlAllQ] 


eW2A2lQl 


cWiai2q2 

W2A22Q2+eW2A22Q2 


(2 . 80) 


Proof:  (i)  Apply  Theorem  2.2  noting  that 

"Vl  0 

V  = 

0  V, 

span  the  left  null  and  row  spaces  of  Aq,  respectively.  (ii)  is  obtained 
by  keeping  in  mind  that  z(t)  =  0. 

There  are  few  interesting  points  to  be  noted.  From  (2.78)  a 
slow  vector  y^  and  a  fast  vector  z^  are  defined  in  terms  of  area  only; 
similarly  for  X2 .  From  (2.80)  the  fast  variables  zpZ2  are  on^y  weakly 
connected  to  each  other  and  since  VT  A^  Q^,  i=l,2  are  nonsingular  these 
connections  can  be  neglected  for  an  0(e)  approximation.  The  fundamental 
difference  between  the  original  (2.71)  and  the  transformed  system  (2.80) 
is  that  (2.80)  no  longer  has  a  continuum  of  equilibrium  points.  Hence, 
each  area  defines  a  local  fast  model  z^  connected  with  0(e)  connections 
to  other  local  models,  whereas  contributions  y^  from  each  area  form  a 
"slow  core"  describing  the  system-wide  dynamics  of  (2.71).  It  is  shown 
in  Chapters  3  and  4  that  this  decomposition  carries  over  to  nonlinear 
weakly  connected  systems. 


35 


CHAPTER  3 

NONEXPLICIT  SINGULAR  PERTURBATIONS  IN  NONLINEAR  SYSTEMS 

3 . 1  Introduc  tion 

An  asymptotic  procedure  for  time-scale  separation  in  nonlinear 
models  is  of  paramount  importance  since  a  nonlinear  analog  of  the  algebraic 
transformation  in  [37,38,41]  is  not  available.  In  this  chapter  we  give  such 
a  procedure  and  demonstrate  its  application  to  classes  of  nonlinear  systems. 
The  coordinate  free  characterization  of  singular  perturbations  [16]  is 
extended  to  nonlinear  systems  of  the  form  ex=h(x, e)  and  it  is  shown  that 
equilibrium  and  conservation  properties  lead  to  a  definition  of  new 
coordinates  in  which  time  scales  are  explicit  (Section  3.2).  In  Section 
3.3  we  study  a  class  of  nonlinear  high  gain  feedback  control  systems  in 
which  the  controls  enter  linearly  through  a  constant  matrix  but  the  open 
loop  system,  the  output  map  and  the  feedback  law  may  all  be  nonlinear.  It 
is  shown  that  these  systems  can  be  studied  through  singular  perturbation 
techniques  after  they  are  transformed  to  the  explicit  form  using  the 
method  of  Section  3.2.  The  last  section,  3.4,  is  devoted  to  interconnected 
systems  whose  isolated  subsystems  posses  equilibrium  and  conservation 
properties.  It  is  shown  that  in  such  systems  weak  connections  give  rise 
to  two-time-scale  behavior.  Separation  of  the  time  scales  defines  a  slow 
"core"  which  describes  the  system-wide  behavior  and  a  set  of  fast  "residues" 
describing  the  local  behavior  of  each  subsystem.  The  decomposition,  known 
for  specific  classes  of  linear  models  such  as  Markov  chains  [42-43]  linearized 
models  of  power  systems  [24-28,44],  electrical  networks  [44,45]  and  economic 


systems  [46]  where  it  appeared  for  the  first  time,  is  established  for  a 
wide  class  of  nonlinear  systems  with  common  features  (i)  equilibrium  and 
conservation  properties  and  (ii)  weak  connections  between  subsystems. 

Using  this  decomposition  we  give  decentralized  stability  cirteria  for  this 
class  of  systems,  analogous  to  those  in  [47-49]. 


3.2  Conservation  and  Equilibrium  Properties  in  Nonlinear  Systems 

The  need  for  coordinate  free  characterization  of  time  scales  in 
nonlinear  systems  is  more  pressing  than  in  linear  systems.  Wide  separation 
of  eigenvalues  provides  some  characterization  in  linear  systems  but  the 
notion  of  modes  is  nonexistent  in  nonlinear  systems.  It  will  be  shown  in 
this  section  that  the  conservation  and  equilibrium  properties  introduced 
in  Chapter  2  for  linear  systems  can  naturally  be  extended  to  nonlinear 
systems  and  that  they  lead  to  a  new  set  of  variables  in  which  the  time 
scales  are  explicit. 

To  motivate  the  discussion  we  re-examine  the  explicit  model  (1.1) 
from  a  different  point  of  view.  Writing  (1.1)  in  the  fast  time  scale 


jZ  =■  ef(y,z,e) 

-  g(y,z,e) 


and  setting  £=*0  we  obtain  the  auxiliary  system 


dz 

dT 


g(y»z.O) 


(3.1) 


(3.2) 


which  has  the  following  two  important  properties. 
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D 


Conservation  Property  A  function  of  the  state 

c( y,z)  =  y  (3.3) 

remains  at  its  initial  value  a(y(0)>  z(0))  =  y(0),  that  is,  it  is  conserved 
during  the  motion  of  (3.2). 

Equilibrium  Property  System  (3.2)  possesses  a  set  of  nonisolated 
(continuum)  equilibrium  points  defined  by 

g(y,z,0)  =  0  .  (3.4) 

The  equilibria  defined  by  (3.4)  are  the  "quasi-steady-states"  to  which  the 
fast  transients  of  (3.1)  converge  as  explained  in  Chapter  1. 

A  generalized  version  of  (3.1)  is  a  system  in  the  form 

e  l?  =  ir  =  h<x’e>  <3-5> 

which  in  t- scale  at  e=0 

=  h(x,0)  (3.6) 

has  equilibrium  and  conservation  properties  analogous  to  the  properties 
of  (3.1).  System  (3.5)  is  studied  in  a  domain  DcRnx[0,eo]  in  which 
function  h  is  assumed  to  be  continuously  differentiable  with  respect  to 
x  and  e. 


Assumption  3.1  System  (3.5)  satifies  the  following  conditions  for 

(*) 

existence  of  manifolds'  S  and  F. 

Equilibrium  Manifold  S  The  set 


S  =  {x|h(x,0)=0,  x£ D  } 


(3.7) 


defines  a  v“dimensional  differentiable  manifold,  v  >  1.  Hence,  there 

exists  continuously  differentiable  function  cp:Rn— R^3,  p=n-v,  rank  ^2=p, 

ox 

V  xCD  such  that 


cp(x)=0  <=>  h(x,0)=0 


(3.8) 


that  is,  in  the  domain  of  interest  D,  every  equilibrium  of  (3.6)  satisfies 

cp(x)=0  and  every  x  satisfying  cp(x)=0  is  an  equilibrium  of  (3.6). 

Dynamic  Manifold  F  There  exists  continuously  differentiable  function 
xo 

ar:Rn— RV  such  that  for  each  x(0)=x  the  p- dimensional  (p=n-y)  manifold 


F  =  l"x|o(x)-a(x  )=0,  rank  ^  =  v] 

xn  U  Ox 


(3.9) 


in  an  invariant  manifold  of  (3.6)  that  is  a  trajectory  originating  in  F^ 
remains  in  F 


cj(x(t))  -  o (Xq)=0*  V  t  >  0  . 


(3.10) 


^Manifolds  are  generalizations  to  Rn  of  objects  such  as  curves 
and  surfaces  in  R  .  More  precisely,  'let  §:Rn-R  be  a  continuously  differenti 

able  function  from  Rn  into  Rm.  Then  if  the  set  M  =  {x|§(x)=0  and  $(x)  has 
rank  m}  is  nonempty,  it  is  an  r-dimensional  manifold,  r=n-m  [50,51]. 
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Moreover,  for  all  x  €D,  manifolds  S  and  F  are  not  tangent  to 

°  xo 

each  other,  that  is,  for  all  x  in  the  intersection  of  S  and  F 


(3.11) 


Theorem  3.2  Under  assumption  3.1,  the  change  of  coordinates 


y=cr(x) 


z=cp(x) 


(3.12) 


transforms  (3.5)  into  the  separated  explicit  model  (1.1)  with 

dz  e=0 

nonsingular  and  'z(t)  =  0. 

Proof:  Differentiating  (3.10)  and  using  (3.6)  we  have 


ff  h<*-°>  *  0  • 


(3.13) 


Differentiating  y=cr(x)  and  using  the  mean  value  theorem  in  e  for  each 
component  of  h 


*1  =  I  32  wx  =  2S.  & 

dt  e  ax  rHX>e;  ax  ae 


(3.14) 


we  see  that  y  is  the  slow  variable  of  (l.l).  Using  the  inverse  transforma¬ 
tion  x*y(y,z)  of  (3.12)  which  exists  because  of  (3.11),  and  differentiating 
z=acp(x)  we  obtain 


dt  “  fx  h(x’e>  =  h(y(y,z),e)  =  g(y,z,e)  . 


(3.15) 


0 


m 


i 


We  show  that  (dg/3z)  is  nonsingular  by  contradiction.  Assuming  that 

e=° 

it  is  singular  the  equilibrium  manifold  of  (3. 14)- (3. 15)  has  dimension 
greater  than  v  which  is  contradiction  because  (3.12)  is  a  nonsingular 
transformation.  Finally  from  x€  S  <=>  x=y(y>0)  it  follows  that  h(Y(y>0),0)=0 
and 


g(y,0,0)=0  (3.16) 

implying  that  "z(t)  =  0. 

The  intuitive  idea  behind  this  theorem  is  illustrated  by  Fig.  3.1. 
If  the  equilibrium  manifold  S  is  attractive,  the  trajectories  of  (3.6) 
which  are  confined  to  some  F  due  to  (3.10),  converge  to  S  and  when  t— 00 
they  terminate  at  the  intersection  of  F  and  S.  Instead  the  trajectories 
of  (3.5)  rapidly  approach  S  staying  in  a  boundary  layer  close  to  F  and 
then  slowly  continue  their  motion  remaining  close  to  S.  Since  the 
trajectories  are  initially  close  to  F  the  quantity  cr(x)  stays  almost 
constant  during  this  interval;  thus  it  qualifies  as  a  predominantly  slow 
variable.  On  the  other  hand,  the  quantity  cp(x),  which  is  large  away  from 
S  where  the  trajectory  starts,  rapidly  diminishes  when  the  trajectory 
approaches  S;  thus  it  qualifies  as  a  predominantly  fast  variable. 


As  an  illustration  we  consider  (3.17) 

d*i 

~  e  -  cp]/x)  +  (XL+X3>  92<x>  "  exi 

dx2  3 
~  =  "2X2  Cf)2(x)  “  eX2 

dx3 

«  9l(x)  +  (Xj+x3)  cp2(x)  -  ex3 


(3.17) 


r . 
i 
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over  R  3  D  =  {(x^x^x^l  >  1,  x2  >  0.5,  x3  >  0.5}  where  cp^x) , cp2  (x) 

are  continuously  differentiable  functions  defined  over  D.  Setting  e=0 
in  (3.17)  we  obtain 


r 


> 


dx. 


-jp-  =  ’  cp1(x)  +  (x1+x3)  cp2(x) 


dx„ 


"37  =  "2X2  *2W 


(3.18) 


dx„ 


-57  ■  <pL(x)  +  (XL+X3)  92(x) 


for  which 


cpL(x)= 


=0 


92(x)=0 


(3.19) 


define  the  equilibrium  manifold  S.  It  is  easily  verified  that  the 
dynamic  manifolds  are  defined  by  cr(x)  =  a(x(0))  where 

o(x)  =  (xl+x3)x2  . 

The  equilibrium  manifold  and  a  dynamic  manifold  of  this  system  are 
shown  in  Fig.  3.1  where  functions  9^>92  were  chosen  as 


(3.20) 


cpj/x)  =  xi_x3 


cp2(x) 


vwL  ■ 


(3.21) 


3.3  High  Gain  Feedback  and  Disturbance  Rejection 

Use  of  high  gain  in  feedback  loops  has  been  known  to  reduce  the 
effects  of  disturbances,  parameter  variations  and  distortions  [19,52]. 
Early  investigations  using  root  locus  techniques  [53]  have  shown  that 


.'--j 


t 


under  high  gain  some  poles  of  the  closed  loop  system  tend  to  infinity,  a 
characteristic  of  singularly  perturbed  systems.  Similar  behavior  is 
exhibited  by  multivariable  systems  [54].  An  extensive  study  of  high  gain 
in  Linear  Time  Invariant  (LTI)  systems  was  undertaken  in  [17,55]  where  it 
was  pointed  out  that  every  high  gain  system  is  a  singularly  perturbed  one 
and  vice  versa.  In  this  section  we  show  that  the  relation  between  high 
gain  feedback  and  singular  perturbations  extends  beyond  the  class  of  LTI 


systems . 


We  consider  the  system 


x  =  f(x)+Bu 

y  =  g(x) 


(3.22) 


under  the  output  feedback 


u  =  jk(y) 


(3.23) 


and  we  study  the  behavior  of  the  closed  loop  system  when  the  gain  . 

e 

The  state  x  is  a  n-vector  and  the  input  and  output  vectors  both  have 
dimension  m  <  n.  Functions  f,g  are  defined  in  a  domain  DxCRnand  function 
k  is  defined  in  a  domain  D^c:Rm.  All  functions  are  assumed  to  be  differenti¬ 
able  a  sufficient  number  of  times.  Moreover  we  make  the  following  basic 
assumption. 

Assumption  3.3  (a)  Matrix  B  is  full  rank.  Hence,  there  exists  a  vxn 

(v=n-m)  matrix  V,  rank  V=v,  such  that  VB=0. 

“fc 

(b)  There  exists  a  unique  set  point  y  6  Dy  such  that 


k(y  )  -  0 


(3.24) 


C 


) 


(c)  Matrix  gxJ  ,  where  is  the  partial  derivative  of  g,  is 


nonsingular  Y  xg  satisfying  g(x)-y  =  0. 


Fart  b  of  the  assumption  could  be  relaxed  to  allow  isolated  roots 
of  (3.24).  However,  nonisolated  roots  which  imply  dead  zones  are  excluded. 
When  the  output  is  linear  in  x,  y=Cx,  Assumption  3.3(c)  requires  to  b 

nonsingular  which  is  equivalent  to  the  assumption  of  [17,55]  that  CB  is 
nonsingular.  Indeed  if  Assumption  3.3(c)  holds  [c]  . •  M  must  be 


full  rank  and  (CB)  1  exists;  conversely  if  (CB)  1  exists  the  row  space  of 


C  and  the  left  null  space  of  B  are  disjoint  and  [c]  is  nonsingular. 
Theorem  3.4  Under  Assumption  3.3  the  transformation 


y  =  Vx 
Js 


(*) 


z  =  g(x)-y 


(3.25) 


transforms  the  high  gain  system  (3.22)- (3.23)  into  an  explicit  singularly 

_  \ 
perturbed  form  with  z=0. 

Proof:  Substituting  (3.23)  into  (3.22)  and  rescaling  the  time,  T=t/e  we 
obtain 

~  ■  «*(*)  +  Bk(g(x))  (3.26) 

whose  auxiliary  system  is 


3*  =  Bk(g(x)) 


(3.27) 


* 

We  temporarily  change  notation  letting  v  instead  of  y  denote 
the  slow  variable.  The  latter  has  been  reserved  tor  its  more  traditional 
use  as  an  output  variable. 


Since  by  Assumption  3.3(a) 


the  relation 


=  VBk(g(x))  =  0 


(3. 


ct(x)  =Vx,  o(x)=ef  (x(0))  (3. 

defines  the  family  of  m-dimensional  dynamic  manifolds.  By  Assumption 
3.3(b)  all  x£  satisfying 

g(x)=y*  (3. 


are  equilibria  of  (3.27).  Assumption  3.3(c)  implies  that  (3.30)  defines 
a  (n-m) -dimensional  equilibrium  manifold  which  is  transversal  to  the 
dynamic  manifolds.  The  transformation  is  simply  an  application  of 
Theorem  3.2. 

In  the  new  coordinates  (3.22)- (3.23)  becomes 


dz 

dt 


dys 

■dT  "  v£(y(Vz)> 

=  egyf(Y(ys,z)  +  gyBk(z+y*) 


(3.: 


where  x=Y(yy,z)  is  the  inverse  transformation  to  (3.25)  which  exists 
because  of  Assumption  3.3(c). 

Corollary  3.5  If  (i)  Assumption  3.3  is  satisfied  and  (ii)  the  boundary 
layer  system 


Bk(z+y*) 


(3.: 


is  asymptotically  stable,  the  response  of  (3.22)-(3.23)  is  approximated 
over  a  bounded  interval  [0,  T] ,  by 


x(t)  =  v(ys(t),  z(j)  +  0(e) 
y(t)  =  z(j)  +  y*  +  0(e) 

where  y  (t)  satisfies  the  reduced  system 
s 


*1 

(3.33) 


dy 

~  =  vf(Y(ys,o))  •  (3.34)  : 

Proof:  Follows  immediately  from  Theorem  3.4  and  standard  singular  J 

M 

perturbation  results. 

From  (3.33)  the  output  differs  from  the  set  point  y  by  the 

predominantly  fast  variable  "z+0(e).  Hence,  any  disturbance  that  can  be  k 

* 

modelled  as  initial  condition  will  appear  in  the  output  only  over  a  short 
initial  interval.  Note,  however,  that  with  the  assumptions  in  Corollary 
3.5  the  approximation  (3.33)  is  valid  only  over  a  bounded  time  interval 

« 

[0,  T] .  Under  stronger  conditions,  which  essentially  amount  to 
stability  requirements  on  the  slow  system  (3.34),  the  approximation  is 
valid  over  unbounded  intervals  [36]  and  disturbance  rejection  is  indeed 
achieved. 

Note  also  that  the  dynamic  manifold  defined  in  (3.29)  is  linear 
because  we  assumed  that  the  input  enters  linearly  through  a  constant 
matrix.  If  in  addition  the  output  is  linear  in  x,  the  equilibrium  mani¬ 
fold  is  linear  too  and  (3.25)  is  a  linear  transformation.  In  this  case 
the  boundary  layer  system  (3.32)  depends  on  'z  only  and  application  of 


Theorem  3.4  and  Corollary  3.5  is  greatly  facilitated  since  transformation 
(3.25)  can  be  inverted  explicitly  and  stability  of  (3.32)  can  be  checked 
more  easily.  The  following  example  illustrates  the  discussion  above. 
Example  3.6  The  output  of 

*1  =  fl(x)+u 
*2  =  f2^+u 
-  2 

y  =  xL+x2 

is  fed  back  through  the  high  gain  law 

u  =  ^k(y)  =  -~-(-y-y3)  .  (3.3 

This  system  satisfies  Assumption  3.3(a),  3.3(b)  and 

V  -  [1  -1],  y*=0  .  (3.3J 


x  = 


(3.3. 

(3.3> 


It  also  satisfies  Assumption  3.3(c)  if  we  restrict  the  region  of  validity 
to  x^  >  0.6.  Transformation  (3.25)  becomes 


with  inverse 


yg  =  xrx 
2 

z  =  x^+x 


2 

2 


-1  +  V  144  (y  +z) 

x  =  _ 2 _ 

1  2 

* 

-1  +  J  l+4(y+z) 


(3.3' 


(3.4 


and  the  transformed  system  is 


dy. 


“jf  -  tf1(x)-f2(x)| 


x=v(ya»z) 


dz 


6  It  =  eI2xifl(x)+f2(x)] 


x=v(ys,z) 


(3.41) 

[  +  yT+4(y  +z)'  ]  (z+z3)  . 


Since  the  boundary  layer  system 


£7  =  -  /^y  +z>  (z+z3) 


(3.42) 


is  asymptotically  stable  in  the  region  of  validity,  approximation  (3.33) 
holds  for  large  enough  gain. 

Transformation  (3.25)  can  also  be  used  to  analyze  disturbance 
rejection  when  the  disturbance  is  modelled  as  an  input.  Consider 


x  =  f(x,w(t))  +  Bu 

y  ■  g(x) 


(3.43) 


under  the  output  feedback 

«  =  “k(y>  (3.44) 

whose  difference  from  (3.22)-(3.23)  is  that  the  disturbance  w(t)  appears 
as  input  to  the  system.  Substituting  (3.44)  into  (3.43)  we  obtain 

dx 

e  =■  ef(x,w(t))+  Bk<g(x))  (3.45) 

which,  under  Assumption  3.3,  can  be  transformed  to 
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Vf(y(y  ,z),  w(t))^  f(t,y  ,z) 
s  s 


e  ff  =  egxf(Y(ys>z)>w(t))  +  8xBk(z+y*)  =  g(t,yg,z,e) 


(3.46) 


using  (3.25).  In  the  following  corollary  we  make  use  of  a  theorem  in  [36] 
which,  for  convenience  is  reproduced  in  Appendix  I. 

A 

Corollary  3.7  Assume  that  w(t)  is  such  that  f,g,  satisfy  all  the 

conditions  of  the  theorem  in  [36].  Then  under  Assumption  3.3  and  for  high 

1  * 
enough  gain  —  the  output  of  (3.43)  remains  0(e)  close  to  the  set  point  y 

for  all  t£  [0,  »)  . 

Proof :  Since  the  boundary  layer  system  of  (3.46) 

ic 

=  gxBk(z+y  )  (3.47) 

has  z=0  as  its  unique  equilibrium  and  y(0)=y  implies  z(0)=y(0)-y  =0, 
the  fast  part  of  z  is  ”z(t)=0  V  t£[0,<»).  Furthermore,  "z(t)=0.  Hence, 

y(t)  =  z(t)  +  y*  =  y*  +  0(e)  (3.48) 

by  the  theorem  in  [36], 

The  practical  use  of  the  above  corollary  may  be  limited  since  it 
assumes  that  w(t)  is  known  so  that  the  assumptions  in  [36]  can  be  checked. 

A  much  more  desirable  result  would  be  to  establish  (3.48)  for  a  class  of 
inputs  w(t).  Such  a  result  should  draw  upon  the  specific  way  in  which  the 
disturbance  enters  into  the  problem. 


An  appealing  approach  to  large  scale  system  analysis  and  design 


is  to  view  the  system  as  a  collection  of  dynamic  subsystems  interacting 
through  static  interconnections.  The  object  then  is  to  analyze  the 
stability  of  [47-49]  or  design  control  laws  for  the  system  [56]  in  a 
decentralized  fashion,  that  is,  by  testing  the  stability  of  the  subsystems 
or  designing  feedback  control  using  only  the  subsystem  states  or  outputs. 
This  approach  is  based  on  the  premise  that  the  connections  between  sub¬ 
systems  are  "weak"  compared  to  the  internal  connections;  hence,  qualitative 
properties  and  control  design  can  be  performed  on  the  subsystem  level. 

In  this  section  we  show  that  when  subsystems  have  equilibrium 
and  conservation  properties,  weak  connections  give  rise  to  two- time-scale 
behavior.  Subsystems  are  weakly  coupled  in  the  fast  time  scale  but  are 
strongly  coupled  in  the  slow  one.  We  start  by  showing  that  such  cases 
arise  naturally  in  high  gain  decentralized  output  feedback.  Consider  the 
interconnected  system 

*i  =  ^i<xi)  +  Vx)  +  Biui 

yi  =  8i(xi)  i-l,...,i  (3.49) 

r  T  T,T 

x  =  [xL  ...  x^] 

with  decentralized  output  feedback 


We  use  the  word  "connection"  to  mean  physical,  static  interaction 
and  the  word  "coupling"  to  imply  dynamic  interaction. 
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ui  =  \ 


i=l, . . . ,1 


(3.50) 


■where  B^,  g^,  k^  satisfy  Assumption  3.3,  i-l,...,X.  Substituting  (3.50) 
into  (3.49)  and  rescaling  the  time  t  =  ~  we  obtain 


"dr"  =  Bi  ki<8i(xi))  +  £(fi(xi^+hi(x^  i=l,...,J& 


(3.51) 


whose  sub sub terns 


—  =  B.  k. (g.  (x.)) 

dT  l  i'°iv 


i~l, * . . , fa 


(3.52) 


have  equilibrium  and  conservation  properties  as  argued  in  section  3.3. 


Generalizing  (3.51)  we  consider  weakly  interconnected  systems 


of  the  form 


-jj-  =  f1(xi,e)  +  egi(x,e) 


i-1 ,  •  • . ,  X 


(3.53) 


“i 

where  f,  is  defined  on  a  domain  D.x[0,e.]cR  x  R,  g  is  defined  on 
1  i  i .  x 

m  X  T  T  * 

Dx[0,e]cR  x  R,  x=[x  ...  x  ]  and  n  =  E  n. .  Functions  f  and  g.  are 

1  /  i=l  ^  it 

assumed  to  be  sufficiently  smooth.  In  (3.53)  f^(x.,e)  represents  the  ith 
isolated  subsystem  [47]  whereas,  eg^(x,e)  represents  interconnections  with 
other  subsystems.  Although  the  dependence  of  f^  on  e  may  seem  superfluous, 
it  is  sometimes  needed  to  assure  the  existence  of  equilibria  of  the  isolated 
subsystems;  such  a  case  is  discussed  in  the  next  chapter  in  relation  to 
dynamic  networks.  Concerning  the  isolated  subsystems  we  make  the  following 
assumption. 

Assumption  3.8  Every  isolated  subsystem 


i=l , . . . ,  JL 


fi(xi,0) 


(3.54) 
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has  equilibrium  and  conservation  properties.  That  is 
(i)  the  set 


Si  =  Cxi|fi(xi,0)=0} 


(3.55) 


is  a  v^-dimensional  equilibrium  manifold  of  (3.54),  0  <  <  n^  ;  hence, 

ni  Pi 

there  exists  smooth  function  cp^:R  such  that 

cpi(Xi)=°  fi(xi,0)=0  (3.56) 

ni  vi 

(ii)  there  exists  function  a^:R  "*  R  such  that 

Ft  =  {xi|ai(xi)  =  ci(xi(0))}  (3.57) 

is  a  family  of  invariant  manifolds  of  (3.54)  parametrized  on  o^(x^(0)). 
Moreover,  and  are  nontangent,  i.e. 


rank 


^ix 

aix 


n 


V  x.ec. 

i  l 


(3.58) 


where  co.  .a.  are  the  Jacobian  matrices  of  cp,  ,c,  . 

ix  Ti  i 

Corollary  3.9  Under  Assumption  3.8  the  interconnected  system  (3.53)  is 
a  two -time-scale  system  and  the  transformation 


zi  = 


i=l, . . . ,1  (3.59) 


A 

transforms  (3.53)  into  an  explicit  model  with  v  =  E  v.  predominantly  slow 

A  i=l  1 

variables  y  and  p  =  £  p  predominantly  fast  variables  z  for  which  z=0. 

i=l  1 


53 


Proof:  Isolated  subsystems  (3.54)  form  the  auxiliary  system  of  (3.53) 
obtained  by  setting  e=0. 


we  obtain  equilibrium  manifold  cp(x)=0,  and  dynamic  manifolds  (j(x(t))=o‘(x(0)) 

of  the  interconnected  system  satisfying  Assumption  3.1.  The  conclusion  is 

then  an  application  of  Theorem  3.2. 

None  that  the  dimension  of  manifold  Si  of  each  subsystem  is 

allowed  to  take  extreme  values  0  and  n^.  When  v^=0  (3.54)  has  at  the  most 

ni 

isolated  equilibria  and  its  dynamic  manifold  is  the  whole  space  R  ;  when 

ni 

v^=n^  (3.54)  is  made  of  n^  integrators  and  its  equilibrium  manifold  is  R 
Note  also  that  transformation  (3.59)  is  block  diagonal  in  the  sense  that 
y^  and  are  defined  in  terms  of  subsystem  state  x^  only. 

Noting  that 


*i  ■  YiOvV 

x  =  y(y,z)  =  [y(y1,z1)  .  .  .  y^y £,z^)] 


(3.61) 


is  the  inverse  transformation  of  (3.59)  which  exists  due  to  (3.58)  and 
rescaling  the  time  t=£T  the  transformed  model  is 
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'ix 


FdfiCYiCyiiZi)’6)  ”1  a 

I  - — -  +  gi(Y(y»z)»e)J  "  Fi(y,z,e) 


dz 

=  cpix  [fi(Yi(yi»zi),0+egi(Y(y»z)»e)l 


4 


G. (y. ,z  ,e)+eH  (y,z,e) 

X  1  1  X 


(3.62) 


where  aix>9^x  are  the  partial  derivatives  of  CT^,cp^  with  respect  to  x^ 
According  to  Corollary  3.9  the  quasi-steady-state  is  "zK),  and  the  slow 
model  is 

d^i 

=  Fi(y,0,0)  i-l,...,l  •  (3.63) 

Rescaling  back  to  t  and  setting  e=0  in  the  second  equation  of  (3.62)  we 
obtain  the  fast  model 

dz. 

-j—  "  Gi(yi,zjL,o)  i“l» . . .  ,X  (3.64) 

where  y^^  appears  as  a  parameter. 

Note  that  is  a  function  of  the  whole  y  vector  whereas,  is 

a  function  of  z^  only.  Hence,  separation  of  time  scales  has  resulted  in 

a  decomposition  in  which  parts  from  every  subsystem  are  put  together  to 

form  a  slow  core  (y-variables)  while  the  rest  of  each  subsystem  forms  a 

fast  residue  (z^-variables) .  The  slow  core  describes  the  system-wide 

dynamics  which  due  to  the  weak  connections  between  subsystems  become 

significant  only  in  the  long  run.  The  fast  residues  describe  the  local 

# 

dynamics  which,  due  to  the  strong  connections  within  subsystems,  are 
significant  in  the  short  run.  If  further,  the  fast  residues  are  asymptoti¬ 
cally  stable  the  z^  variables  reach  quickly  their  quasi-steady-state 


equilibrium  (z^O)  ;  hence,  they  are  weakly  coupled  with  each  other  since 
interaction  from  subsystem  to  subsystem  through  the  weak  connections 
becomes  noticeable  only  in  the  slow  time  scale.  (See  (3.64)).  Figure  3.2 
gives  a  pictorial  view  of  the  discussion.  This  decomposition  is  very 
reminiscent  of  Simon  and  Ando's  reasoning  in  their  classical  1961  paper 
[46] .  We  quote: 

(1)  We  can  somehow  classify  all  the  variables  in  the  economy 
into  a  small  number  of  groups; 

(2)  We  can  study  the  interactions  within  the  groups  as  though 
the  interaction  among  groups  did  not  exist; 

(3)  We  can  define  indices  representing  groups  and  study  the 
interaction  among  these  indices  without  regard  to  the  interactions  within 
each  group. 

Step  (1)  corresponds  in  our  case  to  identifying  subsystems 
connected  to  each  other  through  weak  connections.  Step  (2)  corresponds 
to  our  fast  models  (3.64)  which  are  disconnected;  we  went  one  step  further 
to  remove  the  slow  motion  from  each  subsystem.  Step  (3)  corresponds  to 
the  definition  of  slow  variables  y^  as  "indices"  representing  subsystems 
and  the  study  of  the  system-wide  dynamics  through  the  slow  core  (3.63). 

In  the  next  chapter  where  the  decomposition  is  specialized  to  dynamic 
networks,  the  slow  "indices"  take  the  meaning  of  aggregate  physical 
variables. 

A 

There  is  an  extensive  literature  devoted  to  stability  analysis 
of  interconnected  systems  [47-49  and  references  therein] .  The  general 
plan  followed  is  (i)  to  regard  the  large  scale  system  as  an  interconnection 
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of  isolated  subsystems,  (ii)  to  characterize  stability  properties  of 

isolated  subsystems  through  Lyapunov'  techniques,  (iii)  to  deduce  stability 

properties  of  the  overall  system  from  stability  of  subsystems  and  the 

nature  of  interconnections.  A  basic  assumption  in  [47,48]  is  that  the 

subsystems  have  isolated  equilibria  and  the  results  of  [47,48]  do  not 

directly  apply  when  the  subsystems  have  nonisolated  (i.e.,  a  continuum  of) 

equilibrium  points  such  as  in  (3.53)  with  f^x^jO)  satisfying  Assumption  3.8. 

We  now  show  that  stability  criteria  analogous  to  those  in  [47]  can  be 

derived  based  on  recent  results  [20]  on  stability  of  singularly  perturbed 

systems.  For  convenience,  these  results  are  reproduced  in  Appendix  II. 

We  consider  the  interconnected  system  (3.53)  when  f^(xi,0) 

satisfies  Assumption  3.8.  Using  transformation  (3.59)  we  obtain  the  system 

in  the  explicit  singularly  perturbed  form  (3.62).  We  assume  that 
vi  pi 

y.€D  cR  ,  z.€D  cR  ,  i=l,...,i  and  that  y=0,z=0  is  the  unique 
1  yi  1  zi 

equilibrium  of  (3.62).  Moreover  the  following  assumptions  are  made 
concerning  the  slow  core,  the  fast  residues  and  their  interactions. 

Assumption  3.10 

(i)  The  slow  core  (3.63)  has  a  Lyapunov  function  V:Rv—  R+ such  that  for  all 

y^Dy 

[VyV(y)]T  F(y,0,0)  <  -  aiY2(y),  ^  >  0 


where  y(y)  is  a  scalar  valued  function  of  y  with  y(0)=0  and  Y(y)#0,  y^O. 

(ii)  Every  isolated  fast  residue  (3.64)  has  a  Lyapunov  function 
vi  pi 

W,  (y. ,  z  )  :R  xR  -*  R  such  that  for  all  y.  6D  ,z,€D 
iii  +  i  y^  i  z^ 


[?ZiWi(yi>zi)JT  Gi(y1,zi,0)  <  -  a2t  (zi>  <*21  >  0. 


(iii)  There  exists  Xi  >  0,  i=l,...,jJ  such  that  for  all  y€Dy,  z£Dz  the 
following  hold: 


(a)  [VyiW.(yi,Z.)]T  F.(y,z,e)  <  c.  z  \t  ^(z.) 

(b)  [v  V (y) ]T  [F(y,z,e)-F(y,0,e)]  <3  ¥(y)  (E  \±  ^(z.))L/2 

(c)  [Vz  Wi(yi,zi)]T  [Gi(yi,z.Je)-Gi(yi)zi,0)]  <  e  ku  (J)^(zi) 

(d)  [Vz  Wi(yi,zj.)]T  Hi(y,z,e)  <  k2  E  X± 

i  i=l 

Constants  c^,<S,  k^,k2  are  all  assumed  to  be  nonhegative  numbers. 

Theorem  3.11  If  Assumption  3.10  is  satisfied  and  e  is  sufficiently  small 
the  equilibrium  (y=0,z=0)  of  (3.62)  is  asymptotically  stable. 

Proof:  Let 


J i 

W(y,z)  =  E  \±  W.(y.,z.) 


(3.6: 


be  a  tentative  Lyapunov  function  for  the  boundary  layer  system,  formed  by 
the  z^-sys  terns,  i  i,  in  (3.64).  Then, 


[V  W(y, z) ]  G(y,z,0)  =  E  ?z  ^t(y  tz  )  GjL(yi,zi,0) 

i=l  i 


(3. 6i 


<  -  2  \i  a2i  ^i(zi)  <  "  a2 

i=l 


where 

(}>2(z)  =  Z  ^  C^(z.) 
i 


^“roin  cy2i 
i 


(3.6 


Hence,  V(y),W(y,z)  satisfy  condition  (I),  (II)  of  [20].  Condition  Ilia 
of  [20]  is  also  satisfied  since 
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[V  W(y,z)]T  F(y, z, e)  =  £  [7  W._(yi’zi)]T  Fi(y,z,e) 

y  i  *i 

<EXici  C^2  (z)  =  (X  Xici)  4*2  (z) 


(3.68) 


whereas  condition  Illb  is  identical  to  (iiib)  in  Assumption  3.10. 
Finally,  letting 


G±(y,z,e)  =  Gi(yi,zi,e)  +  e  Hi(y,z,e) 


(3.69) 


we  have 


T 

i  r 


[7  W(y,z)]  [G(y, z, e)  -  G(y,z,0)]  = 


£  [V  W  (y  ,z  )]T  [G. (y  ,z  ,e)  -  G.(y.,7  ,0)  +  e  H  (y,z,e)]  <  (3.70) 

.  z .  l  i  J.  rir  iir  j. 

i  i 

t  E  ^(z.)  +  e  k2  C|)2(z)  <  e  C|)2  (z) 


where 


k  . 

K.  =  max  •— 1  +  k 
1  i  *i 


(3.71) 


and  we  see  that  condition  IIIc  of  [20]  is  also  satisfied.  The  conclusion 
follows  directly  from  [20] . 

A  similar  procedure  cannot  be  applied  to  the  original  model  (3.53) 
because  the  isolated  systems  possess  a  continuum  of  equilibrium  points.  The 
basic  difference  between  (3.53)  and  (3.62)  is  that  in  (3.62)  all  the  slow 
dynamics  giving  rise  to  equilibrium  manifolds  have  been  relegated  to  the 
slow  core;  consequently,  the  fast  residues  no  longer  have  a  continuum  of 
equilibrium  points.  The  stability  criteria  become  easy  to  apply  when  the 
transformed  models  are  structured,  such  as  in  dynamic  networks,  the  subject 


of  the  next  chapter. 


CHAPTER  4 


REDUCED  ORDER  MODELING  OF  DYNAMIC  NETWORKS 

4.1  Introduction 

Much  of  the  equilibrium-conservation  reasoning  presented  in 
Chapters  2  and  3  was  inspired  by  the  study  of  time  scales  in  power  systems 
and  other  weakly  connected  networks.  In  this  class  of  systems,  which  is 
the  subject  of  the  present  chapter,  states  and  connections  have  physical 
meaning  and  separation  of  time  scales  is  related  to  physical  laws  such  as 
conservation  of  mass,  charge,  momentum,  etc. 

Section  4.2  discusses  weakly  connected  networks  with  linear 
storage  but  nonlinear  interconnection  elements.  The  main  result  is  that 
the  transformation  that  brings  the  system  into  the  explicit  singularly 
perturbed  form  is  linear.  In  the  new  coordinates  a  slow  core  describes 
the  system-wide  behavior  while  fast  residues  describe  the  local  behavior 
of  the  network.  The  slow  core  turns  out  to  be  another  "aggregate"  network 
whose  states  and  connections  are  related  in  an  intuitively  appealing  way 
to  the  states  and  connections  of  the  original.  These  results  are  then 
specialized  to  power  systems  and  a  five-machine  example  illustrates  the 
reduction  procedure  (section  4.3);  simulation  results  are  also  shown. 

In  section  4.4  we  clarify  the  relation  between  coherency  and  localizability 
[24-25]  for  LTI  systems. 
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4.2  Time  Scales  in  Nonlinear  Dynamic  Networks 

The  time-scale  separation  methodology  developed  in  the  previous 
chapter  will  now  be  applied  to  nonlinear  dynamic  networks,  a  class  of 
large  systems  whose  structure  facilitates  the  derivation  of  reduced  models 
with  physical  meaning.  The  dynamic  networks  we  consider  are  systems 
comprised  of  storage  elements,  capable  of  storing  some  physical  quantity 
and  interconnection  elements  capable  of  transporting  this  quantity  without 
delay.  Examples  of  dynamic  networks  include  power  systems,  where  angular 
momentum  stored  in  the  generators  is  transported  through  transmission  lines, 
R-C  networks  where  charge  in  the  capacitors  is  transported  through  resistors, 
mass-spring  systems,  etc.  The  dynamic  networks  considered  first  have 
storage  elements  with  linear  characteristics  but  their  interconnections 
may  be  nonlinear.  Extension  to  networks  with  nonlinear  storage  elements 
is  indicated  later  in  the  section.  The  rates  of  flow  in  the  interconnections 
are  assumed  to  be  continuously  differentiable  functions  of  the  potential 
differences  across  the  interconnections  satisfying 


WvV  ■  -  £ki<W 


(4.1) 


This  assumption  is  equivalent  to  saying  that  there  are  neither  sources 
nor  sinks  along  the  interconnection.  The  dynamics  of  these  systems  are 
then  modeled  by  either  the  system  of  first  order  equations 


-•  E  I 

i  k€  k. 


fik(xi"V  "  V 


(4.2) 


or  the  system  of  second  order  equations 


••  _  ^  r 

x<  7T  E 

1  i  k 


S  fik(W  - 


h] 


(4.3) 


where  x^m  the  potential  and  inertia  of  the  ith  storage  element,  1^ 
the  net  injection  at  the  ith  element,  fikC)  the  characteristic  of  the 
interconnection  between  elements  i  and  j  and  k^  the  set  of  elements  to 
which  i  is  connected.  In  the  remainder  of  this  section  we  deal  with 
dynamic  networks  in  the  form  (4.2).  In  the  next  section  the  results  will 
be  applied  to  power  systems  whose  equations  are  in  the  form  (4.3). 

A  dynamic  network  is  said  to  be  weakly  connected  if  some 
interconnections  can  be  expressed  as  multiples  of  a  small  parameter  e.  The 
model  of  a  weakly  connected  network  is  then 


dx. 


dr 


k€K, 


fik(VXk> 


+  e 


8ij (xi' 


V 


I^(e)] 


(4.4 


where  K^J.^  are  index  sets  representing  nodes  connected  to  element  i. 

Constant  1^,  which  depends  on  and  is  differentiable  with  respect  to  e,  is  a 

net  injection  (of  power  or  current)  at  node  i.  Its  dependence  on  e  will  be 

discussed  later.  A  fundamental  property  of  weakly  connected  networks  is 

that  neglecting  the  weak  connection  terms  eg^.  results  in  v  isolated  "areas' 

a*>l,...,V’  Area  a  contains  n  connected  nodes  and  its  equation  is  obtained 

a 

by  setting  e=0  in  (4.4) 


i 


fik(xi-xk)  "  Ii<0)J’ 


i €  os  *  •  *  *  j  v*  (4.5) 


When  the  states  in  each  area  are  ordered  consecutively  the  nxv  partition 


matrix  U  is 
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U  ■  diag  (u^,...,u^)  (4.6) 

where  u  is  an  a  -vector  with  all  elements  one. 

a  of 

Assumption  4.1  Each  of  the  v  areas  formed  by  setting  e“0  in  (4.4)  has  an 
equilibrium  state. 

In  the  case  of  power  systems  the  above  assumption  requires  that 
every  area,  when  isolated  from  the  rest  of  the  system,  has  its  own  load 
flow.  This  will  only  be  possible  if  the  area  adjusts  its  net  injections 
I^(e)  so  that  the  power  exchanged  with  other  areas  is  compensated  internally. 
Thus,  the  injections  I^(e)  are  made  functions  of  the  strength  of  the  inter¬ 
area  connections  e.  The  choice  of  the  dependence  of  I ^  on  e  and  its  impact 
on  the  accuracy  of  the  reduced  models  will  be  discussed  later  in  the  section. 
In  each  area  a  we  select  a  reference  node  xr>  r£ & v  * ,  and  form 

the  difference  sir  ■  x®-x®  for  r€<*  and  V  iga  ,  ij*r,  where  x®,x®  are  the 

£ 

values  of  x^«xr  at  an  equilibrium  x  of  the  area  model  (4.5). 

Theorem  4.2  System  (4.5)  has  an  equilibrium  manifold  S  described  by 

cPiOO  ■  xi"xr"sir  =  0  (4.7) 

for  r€<*;  V  i£a,  i^r  and  all  areas  The  dynamic  manifold 

F  for  x(0)=x  is 
x  o 

o 

a  (x)  -  a  (x  )  -  0  ,  (4.8) 

a  a  o 

s 

where 

a  (x)  “  £  b.x,/  E  m  .  (4.9) 

a  iSa  i€» 


v  'Abusing  notation,  we  let  a  be  the  index  of  an  area  as  well  as 
the  set  of  node  indices  in  the  area. 
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Furthermore, 

y  (x)  -  a  (x)  ,  2,  -  ^(x)  (4.10) 

a  a  1  1 

are  v  slow  and  p=n-v  fast  variables  satisfying  Theorem  3.2. 

Proof:  Any  x  satisfying 

*  xi_x£  V  i,k€o,  a=l,  ...»v  (4.11) 

is  an  equilibrium  of  (4.5).  Following  a  path  from  node  i  to  node  r, 
these  relations  can  also  be  written  as 

vv<v*i+i>+<xi+r*i+2>+-  •  •+<v2‘xr-i  )+<vrVVN 

(4.12) 

r€o»  V  16  o,  V  a 


which  is  the  expression  in  (4.7). 

Writing  (4.5)  at  an  equilibrium 


-  ~  [  Z 
i  k€K, 


Lik 


,  e  e. 

<VV 


^<0)1 


l€o,  o=l»*««»v  (4.13) 


and  using  (4.1)  we  obtain 


Z  1,(0)  =  0  ,  o=l,.. .,v  •  (4.14) 

i€  a 

The  last  relation  gives 

2  m1(dxi/dr)  =  0  ,  o“l,...,v  (4.15) 

l€o 

where  (4.1)  was  used  once  more.  The  dynamic  manifold  (4. 8-4. 9)  is  obtained 
by  integrating  and  scaling  (4.15).  Finally  the  transformation  (4.10)  is  an 
application  of  Theorem  3.2. 
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Note  that  although  the  model  is  nonlinear  both  the  equilibrium 
manifold  (4.7)  and  the  dynamic  manifold  (4.8)  are  linear  leading  to  a 
linear  transformation  separating  the  time  scales.  Manifold  S  is  linear 
because  the  right-hand  side  of  (4.5)  is  a  function  of  a  linear  combination 
(the  differences)  of  the  states  as  opposed  to  being  a  function  of  the 
states  individually.  Manifolds  F  are  linear  because  the  conservation 
property  is  linear.  In  the  case  of  RC-circuits  the  conservation  property 
expresses  Kirchhoff's  current  law  (KCL)  and  in  the  case  of  power  systems 
the  conservation  of  angular  momentum.  These  physical  laws  are  linear  even 
when  some  elements  of  the  network  have  nonlinear  characteristics. 

To  rewrite  (4.10)  in  matrix  form  we  define  the  difference 
matrix  G  ■  diag  (G^,...,G^)  where 


-110 
-10  1 


0 

0 


L 


(4.16) 


{.s  a  (n  -l)xn  matrix  with  two  nonzero  elements  per  row.  Ordering  the 
a  a 

states  x  in  the  same  area,  consecutively,  with  the  reference  state  first, 

T  -IT 

denoting  M  *  diag  (m.,...,m  )  and  C  ■  (U  MU)  U  M  the  transformation 

1  cL 

(4.10)  is 


0 

S 


(4.17) 
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where  s  is  a  p-vector  with  components  s^r.  The  inverse  of  (4.17)  is 


x  =  [U  B] 


(4.18) 


-IT  -IT  -1 

where  B  =  M  G  (GM  G  )  .  Recall  from  (4.6)  that  every  row  of  U  has 

one  entry  1  and  the  rest  0,  whence. 


x.  =  y  +bj  (z+s) 

l  a  t. 


V  i^a  >  cr5!*  •  •  • » v 


(4.19) 


where  is  the  i  row  of  B.  After  simple  manipulations  we  obtain  the 


transformed  model 


te  iL  jIj,  8u  1(VV+<VV(2+S)I 

3=1, ...» v 


E  M«>  (4.20) 


57  '  =  £ik  <2rV5lr-8kt)+£lc(zi+Sit) 

i  k6Kt 


S-  t  S  hv  <2k+V>+  V” 

*  r 


*  >7  giJ  [(VV+<W(2+S>1 

~~  2  grj  [(ya-yp)+(br-bj)(z+s)]] 


(4.21) 


In  (4.20),  (4.21),  i  €  or ,  j  €  3  and  b^,  bj  ,  bf  are  the  i,  j,  r  rows  of  B, 
respectively,  r €«*.  Since  from  (4.14)  £  1^0)  =  0,  a=l,...,v 
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2  I  (e)  -  0(e)  ,  a-1, ...,v  (4.22) 

16<y 


and  the  right-hand  side  of  (4.20)  is  0(e). 

In  the  original  state  description  (4.4)  areas  cannot  be  considered 
weakly  coupled  because  over  a  longer  period  their  interaction  through  weak 
connections  becomes  significant.  In  the  transformed  description  the  fast¬ 
time  area  models  in  z  variables  are  weakly  coupled,  and  the  long  term  area 
interaction  is  approximately  described  by  (4.20),  that  is,  the  aggregate 
variables  y  alone.  The  fundamental  difference  between  this  model  and  the 
original  is  that  the  decoupled  z-equations  obtained  by  setting  e-0  in 
(4.21),  no  longer  have  a  continuum  of  equilibrium  points. 

The  definition  of  slow  coherency  as  given  in  [24-28]  is  based 
on  a  modal  decomposition  and  is  not  directly  applicable  to  nonlinear 
systems.  Since  we  have  shown  that  the  two- time- scale  properties  remain 
valid  for  nonlinear  systems,  we  will  use  them  for  the  following  generaliza¬ 
tion  of  the  notion  of  slow  coherency. 

Slow  coherency.  States  x^,xj  (4*4)  are  said  to  be  slow 
coherent  if  x(0)€S  implies  x^(t)-x^(t)  -  const  V  t  >  0.  States  x^,x^ 
are  said  to  be  near  slow  coherent  if  there  exists  a  bounded  function  of 
time  £(t)  such  that 


x(0)€S  =>  xt(t)-Xj(t)  -  const.  +  e£(t),  V  t  >  0.  (4.23) 

An  area  is  slow  coherent  if  any  two  states  in  the  area  are  near  slow- 
coherent.  The  following  theorem  relates  weakly  coupled  areas  and  slow 
coherent  areas  extending  the  corresponding  result  in  [24-28] . 


Theorem  4.3  If  (4.21)  satisfies  Assumption  1.2  and  e  is  sufficiently 
small,  system  (4.4)  has  v  slow-coherent  areas  specified  by  U. 

Proof:  From  Theorem  4.2,  if  i,j€a 


x.-x.  »  (x.-x  )-(x.-x  )  =  z.-z.  +s.  -s  .  . 
ij  '  i  r'  v  j  r'  ij  air  sjr 

If  x(0)6  S  ,  z(0)=cp(x(0))=0,  which  combined  with  (3.16)  and  (1. 

that  z(t)  =  0(e).  Then  (4.23)  follows  from  (4.24). 

Model  (4.20),  (4.21)  is  in  the  explicit  form  (1.1). 
21i(e) 

letting  e-*  0  ,  I  =  lim  — - -  and  using  the  fact  (Theorem  3.2, 

that  ~z  -  0  the  slow  model  is 


dy  /dt  =  -(1/  E  £  gt.((y  -ya)+0>i“b  )s] 
®  j€Jt  ij  “  e  1  J 

0«1, . . . , v 

* 

“  S  Ij  J  Ci~  1».».»\) 

i€  a 


The  fast  model  is 
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(4.24) 
4)  implies 

Hence, 

Theorem  4.2) 

(4.25) 


(4.26) 


i€  or »  i^r  cTlf.fV. 


Slow  model  (4.25)  represents  an  aggregate  dynamic  network  with  storage 
elements 


m 
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net  Injections 
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and  interconnection  characteristics 


<y  -yfl) 

or  0 


iL  ^‘VV  (bi'Vsl- 
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(4.27) 


(4.28) 


(4.29) 


The  aggregate  model  (4.25)  is  decoupled  from  the  local  models 

(4.26) .  Since  the  sums  in  (4.26)  involve  nodes  from  the  same  area,  the  z 
equations  for  two  different  areas  are  decoupled,  that  is,  the  fast  models 

(4.26)  are  local  in  the  sense  that  they  involve  quantities  from  one  area 
only.  Thus,  each  area  uses  its  local  model  and  at  the  same  time  provides 
the  data  and  receives  results  from  the  global  model.  This  multimodeling 
decomposition  helps  in  formulation  of  decentralized  controls  [57]. 

In  (4.4)  and  in  subsequent  derivations  it  was  assumed  that  the 

dependence  of  the  injections  I^(e)  on  e  is  known.  In  a  realistic  situation 

e  has  a  specific  value  and  injections  are  constant.  The  dependence  on  e 

» 

is  an  asymptotic  tool  guaranteeing  that  the  isolated  areas  formed  by 
e-0  have  a  well  defined  equilibrium.  Therefore,  for  any  function  I(E) 


satisfying  y  1,(0)  =  0,  0=1, . .., v  Assumption  4.1  holds  and  the  quantities 

i€a 

in  (4.7)  are  well  defined.  However,  the  freedom  in  choosing  1(e) 

can  be  utilized  to  influence  the  accuracy  of  reduced  models  in  realistic 

systems  where  e  may  not  be  very  small. 

Note  that  for  e  >  0,  the  equilibria  of  (4.25),  (4.26)  are 

generally  different  from  the  equilibria  of  (4.20),  (4.21).  It  has  been 

observed  in  numerical  experiments  that  the  approximation  of  the  time 

response  improves  when  the  equilibrium  of  the  reduced  models  is  closer 

to  the  equilibrium  of  the  original  (4.20),  (4.21).  It  is  desirable  to 

make  the  two  equilibria  as  close  as  possible,  particularly  for  oscillatory 

responses,  and  if  the  reduced  model  (4.25),  (4.26)  is  used  for  stability 

analysis.  The  following  corollaries  provide  guidance  in  this  direction. 

£ 

Corollary  4.4  Let  x  be  the  equilibrium  of  (4.4).  The  equilibria  of 
(4.25),  (4.26)  are  equal  to  the  equilibria  of  (4.20),  (4.21)  if  and  only 
if 

s  =  GxE  .  (4.30) 

Proof:  First  note  that  by  Theorem  3.2  and  Theorem  4.2  z  =  0  is  the 
equilibrium  of  (4.26),  irrespective  of  the  choice  of  s.  From  (4.17) 

E 

the  equilibrium  of  (4.21)  corresponding  to  x  is 

zE  =  GxE-s  (4.31) 

« 

which  is  made  zero  by  (4.30).  Setting  z  =  0  in  (4.20)  to  obtain  (4.25) 

£ 

does  not  alter  the  equilibrium  because  z  =0.  The  choice  (4.30)  is  unique 
because  (4.31)  is  linear  in  s. 


Corollary  4.4  shows  that  there  is  a  unique  s  for  which  the 
equilibria  of  the  exact  and  approximate  systems  are  equal.  Since  by 
definition  s  =  Gx6  and  G  is  a  matrix  of  full  rank,  (4.30)  implies  that 

xe  =  xE  .  (4.32) 


The  next  corollary  gives  a  necessary  condition  on  I^(e)  such  that  (4.32) 
is  satisfied.  Boundary  nodes  are  nodes  to  which  interarea  connections 
are  attached. 

Corollary  4.5  Equation  (4.32)  is  satisfied  only  if 

I^e)  -  Ii(0)  =  e  E  (x®-xE)  (4.33) 

j  6 


that  is,  the  net  injection  at  boundary  nodes  is  adjusted  by  the  interarea 
flow  while  it  is  left  unaltered  at  nonboundary  nodes. 

Proof:  For  the  equilibrium  of  (4.4) 


(Xi-Xfc)  z  £  8ij  (xi"xj)  +  Ii(e) 


(4.34) 


and  for  the  equilibrium  of  (4.5) 


2 

kgK^ 


^(0). 


(4.35) 


Hence,  (4.33)  is  necessary  for  (4.32)  to  be  true. 

In  cases  such  as  water  distribution  networks  where  some  storage 
elements  may  be  nonlinear,  the  reduction  procedure  is  still  applicable 
after  some  modification.  Assuming  that  the  stored  quantity  is  a  strictly 


?????? 
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monotonic  function  of  the  potential,  the  dynamics  of  the  network  are 
described  by  an  equation  analogous  to  (4.2)  in  which  nu  is  now  a  function 
of  x^.  Equation  (4.9)  becomes 


a  (x) 
a 


£  m.  (x.  )x, 
i  l  i 
_ 

i€  ot 


(4.36) 


and  the  dynamic  manifold  is  no  longer  linear.  The  equilibrium  manifold, 
however,  is  still  given  by  (4.7). 


4.3  Power  System  Application 

The  concepts  of  area  aggregation  and  slow  coherency  that  emerged 
from  the  separation  of  time  scales  in  dynamic  networks  (Theorem  4.2) 
originated  as  model  order  reduction  techniques  in  power  systems  [21-23]. 
First,  a  group  of  coherent  generators,  that  is,  generators  that  "swing 
together"  is  identified  and  then  this  group  is  replaced  with  an  equivalent 
generator.  Analytical  studies  of  coherency  [58,59]  and  coherency  based 
aggregation  [24-28]  were  based  on  linearized  versions  of  the  electro¬ 
mechanical  model  of  power  systems.  In  this  section  we  apply  the  reasoning 
of  Section  4.2  to  extend  the  model  simplification  approach  in  [24-28]  to 
the  nonlinear  electromechanical  model,  and  to  more  complex  models  involving 
flux  decay  dynamics  and  voltage  regulator. 

The  well  known  electromechanical  model  [60]  of  multimachine 

» 

power  systems  is 


2H.  5i=  P  .-P  . 

1  1  mi  ei 


i  1,2, ...  ,7] 


(4.37) 
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where  6^  is  the  rotor  angle  of  machine  i,  is  its  inertia  constant, 

p  and  P  are  its  mechanical  input  power  and  electrical  output  power, 
mi  ei 

respectively,  and  the  small  damping  was  ignored  [25] .  In  this  model 
is  assumed  to  be  constant  and  Pg^  is  given  by 


ViVj  Bij  Sia  <6i“6j>  +  vt  Gii 

J# 


(4.38) 


where  v^  is  the  constant  voltage  "behind  the  transient  reactance,"  B„ 
is  the  (ij)-th  entry  of  the  admittance  matrix  reduced  to  the  machine  nodes 
and  represents  the  load  conductance  at  node  i. 

Substituting  (4.38)  into  (4.37)  we  obtain 


•• 


2Hi  ViVj  Bij  Sltl  (6i‘6j)'(Pmi“Vi  Gii)J 


(4.39) 


which  is  in  the  form  of  (4.3)  with  x^-6^. 


fij  (xi~xj)  =  vivj  Bij  sin  (6i“6j) 


Xi  "  Pmi  -  Vi  Gii 


(4.40) 

(4.41) 


Multimachine  power  systems  are  often  comprised  of  groups  of  tightly 
connected  machines  with  weak  connections  joining  the  groups.  Assuming 
that  weak  connections  are  known  system  (4.39)  takes  the  form  (4.4)  for 
which  Theorem  4.2  gives  equilibrium  and  dynamic  manifolds  and  defines  slow 
and  fast  variables.  Note,  however,  that  since  damping  was  neglected  the 
response  of  (4.39)  is  purely  oscillatory  and  the  separation  of  time  scales 


is  understood  in  the  sense  of  separating  low  frequency  from  high  frequency 
oscillations  [61].  Using  trans forma cion  (4.10)  we  can  arrive  at  expressions 
similar  to  (4.20)- (4.21)  from  which  the  slow  core  and  the  fast  residues  are 
defined  as  in  (4. 25)- (4. 26) .  The  reduction  procedure  in  power  systems  and 
some  physical  interpretations  of  the  reduced  models  are  illustrated  through 
the  following  five-machine  example. 

In  the  power  system  of  Fig.  4.1,  =  0.5,  v^_^,  i=l,...,5  and 

B12=B23=B45=1’  B34=B25”B15=B14=0‘ The  net  injections  ^  and  the  resulting 
steady-state  angles  (in  radians)  are  given  in  columns  1  and  2  of  Table  4.1. 

Table  4.1.  Bus  angles  for  five-machine  power  example 


1 

2 

3 

IL  =  -  0.28 

6i-° 

6!  =  ° 

I2  =  -  0.077 

62  =  0.171 

62  =  0.215 

I3  =  0.186 

63  =  0.391 

63  =  0.458 

I.  =  0.362 

4 

6.  =  0.723 

4 

6,  =  1.042 

4 

I5  =  -  0.191 

55  =  0.456 

65  =  0.730 

Note  that  since  admittances  B^,  B15’  B14  are  mucdl  sna^er  than  the 

rest,  the  system  is  divided  into  two  weakly  connected  areas  a  =  [1,2,3], 

3  =  [4,5].  Suppose  now  that  line  is  tripped  and  we  want  to  simulate  the 
resulting  oscillations  using  reduced  models  (4.25)- (4.26) .  The  post  fault 
load  flow  (shown  in  column  3  of  Table >4.1)  gives 
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s  ^2  “  ^^“§2  =  "  0.215 

8 32  -  63-6*  -  0.243  (4.42) 

8 45  *  64-65  =  0.312. 


Defining  new  variables 


5 ^+63+63 


64+65 


(4.43) 


Z1  =  6r62“S12 


63-62“S32 


Z3  *  V65“S45 


and  letting  e— 0  we  obtain  the  slow  model 
•• 

”  0.033  sin  (y^-j^-O.OSS)  -  0.033  sin  (y^-^+0. 147) 

-  0.033  sin  (y^+0. 078) -0.057 

(4.44) 

•  • 

y2  ■  -  0.05  sin  (y2-y1“0.078)  -  0.05  sin  (y  -yj-0. 147) 

-  0.05  sin  (y2-y1+0.068)+0.086 

and  the  fast  model 


«• 

z1  *  -  2  sin  (z^O.215)  -  sin  (?2+0.243)  -  0.185 
•• 

z2  -  -  2  sin  (.z., +0.243)  -  sin  (Zj-0.215)  +  0.269 


(4.45) 


Z3  ■  -  2  sin  (^3+0.312)  +  0.614  .  (4.46) 

» 

Note  that  (4.44)  is  decoupled  from  (4.45),  (4.46).  The  aggregate  model 

(4.44)  represents  the  oscillations  of  the  aggregate  angles  y^,y2  a8®inst 
each  other,  whereas,  the  local  models  (4.45),  (4.46)  represent  the 
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intermachine  oscillations  in  areas  a,g,  respectively.  Figures  4. 2-4. 5 
show  simulation  curves  with  initial  conditions  equal  to  the  prefault 
equilibrium.  Figures  4.2,  4.3  show  exact  (solid  lines)  and  approximate 
(dotted  lines)  responses  of  angles  6^,6^,  whereas.  Figures  4.4,  4.5 
show  exact  (dotted  lines)  and  approximate  (solid  lines)  responses  of  the 
transformed  variables  y^, z^.  Note  that  generator  angles  are  mixed  variables, 
whereas,  y^  is  predominantly  slow  and  is  predominantly  fast. 

We  now  turn  to  more  complex  models  of  power  systems  and  again 
investigate  the  effect  of  weak  connections  on  the  time  scale  behavior  of 
the  system.  The  model  we  employ  is  basically  the  one  in  [62]  with  a  slight 
simplification;  we  do  not  include  the  fictitious  quadrature  axis  coil  g 
which  is  meant  to  model  eddy  currents  in  the  rotor.  With  this  simplification 
the  model  is 


377(w±-l) 


2H1»1 


•  / 

eqi 


“fi 

5fdi 

Si 

Si 


—  -  e'  i  ,  -  D.  (a) -1) 
wt  qi  qi  iv  i  ' 

T{'oi  teqi  "  (xdi_x  di*  idi+Efdi1 
(_B» +  ^  E«i) 

4;  [_(ke  -  sE  <Efdi»  Efdi +  v 

-  2  e^j  cos  (Sj^-Sj) 

j  Bij  eqj  sin  (6i"6j) 


(4.47a) 

(4.47b) 

(4.47c) 

(4.47d) 

(4.47e) 

(4.48a) 

(4.48b) 


C 


Exact  (dotted  line)  and’  approximate  (solid  line) 
responses  of  y, . 


Exact  (dotted  line) and  approximate  (solid  line) 
responses  of  z„. 


*ti  v  'di  'qi  *  'di  di  '  ~  qi  ’  'qi  ~qi 


where  state  e^  is  proportional  to  the  field  flux,  equations  (4.47d-f) 
model  the  voltage  regulator-exciter  system  and  (4.48),  (4.49)  give  the 
interaction  of  the  generators  through  the  transmission  network. 

If  the  power  system  is  made  of  v  weakly  connected  areas,  (4.48a), 
(4.48b)  are  written  as 

*41  ‘  *  ,1  Blj  eqj  COS  <W  *  Blk  *qk  “s  <W  (4-50»> 

jea  or 


Si  =  =  Bij  'qj  Sl"  +  «5  “lk  eqk  Sl"  <W 

H  j€o  J  J  J  Ms 


(4.50b) 


where  i£  a ,  q»=1,  ...,v.  We  now  make  the  important  observation  that  (i) 
machines  interact  solely  through  currents  i^,i^  and  (ii)  i^, i^  are  functions 
of  the  differences  6^—6^  of  angles,  not  of  angles  individually.  Hence, 
letting  e=0  in  (4.50)  and  setting  the  right-hand  side  of  (4.47)  equal  to 
zero  we  see  that  if  an  assumption  analogous  to  4.1  is  met,  points  satisfying 


6i"6r”Sir  -  0 
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■  iHi  .■  pii^VU".  HJ 

\  *1  .  •  -  A  _  ■  - 

-  *  I 

m  . 

be  if—- 


Ni  r 


where  subscript  e  denotes  equilibrium  and  r€ a  »  a »  i^r,  <**l,...,v  are 
equilibrium  points  of  (4.47).  That  is,  (4.47)  at  e*0  has  a  v-dimensional 
equilibrium  manifold  described  by  (4.51).  An  argument  similar  to  the  one 
in  Theorem  4.2  shows  that 

Dt 

c(6,u>)  -  Z  377  61  +  Z  2HiUJi  -  a(6(0),w<0))  (4.52) 

i€a  i-6  a 

for  a*l,..., v  defines  the  family  of  dynamic  manifolds  of  (4.47).  In  all 

realistic  cases  D^«377  so  that  we  can  ignore  the  first  term  in  (4.52). 

Note  that  although  we  started  with  a  higher-dimensional  model  of  the  power 

system  the  equilibrium  and  conservation  relations  (4.51),  (4.52)  involving 

6,(u  variables  are  identical  to  the  ones  obtained  by  working  with  the 

electromechanical  model.  Consequently,  we  obtain  an  electromechanical 

slow  model  involving  the  aggregate  variables  6  *  (  £  2H.  6.)/  £  2H  , 

*  i6c*  i€cr 

to  “6  . 

a  a 

4.4  Coherency  and  Lo^a Usability 

As  shown  in  Section  4,2  weak  connections  in  a  dynamic  network 
give  rise  to  slow-coherent  groups  of  states  which  are  described  by  local 
models.  We  now  restrict  ourselves  to  linear  time  invariant  systems  and 
investigate  the  relation  between  coherency  and  localizability  when  weak 
connections  are  not  present.  This  discussion  clarifies  the  presentation 
in  [24,25]  where  the  two  notions  werp  essentially  treated  as  equivalent. 

For  the  sake  of  completeness  we  repeat  here  the  definition  of  coherency 
of  [24,25]. 


Definition  4.6  Let  x- Ax  be  a  LTI  system,  9  be  a  subspectrum  of  A  and  V 
■  c 

be  the  corresponding  eigenspace.  Then  states  x^  and  x^  are  said  to  be 

coherent  with  respect  to  a  if  and  only  if,  x(0)€V  implies  that 

o 


Xt(t)  -  Xj (t) 


V  t  >  0  . 


(4.53) 


A  group  of  states  is  said  to  be  a  o-coherent  group  if  any  two  states  from 
the  group  are  coherent  with  respect  to  9. 

Suppose  now  that n  states  of  a  system  are  considered  to  belong 
to  a  group  a-  The  criterion  for  such  grouping  can  be  geographic  proximity, 
accessibility  to  remote  sensing  or  similar.  In  an  attempt  to  describe  the 
"local"  behavior  in  the  group  we  use  n^-1  differences  x^-x^,  where  i,k 
belong  to  the  a  set  of  indices.  Typically,  we  fix  index  k  as  the  local 
reference  and  take  all  i#k  in  the  group  <*  to  form  the  differences.  We  then 
investigate  under  what  conditions  the  local  variables 


z  ■  [G  J  0]  x  =  G 


(4.54) 


are  independent  of  the  rest  of  the  system. 
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Definition  4.7  Group  a  is  said  to  be  localizable  if  there  is  an  A^ 

(*) 

such  that 

z  48  A  z  (4.56) 

a 

Lemma  4.8  A  group  of  states  is  localizable  if  and  only  if  it  is  a 
a-coherent  group  and 


n  *  n  -  n  +1  (4.57) 

•a  a 

where  n  is  the  number  of  modes  in  cr. 

<T 

Proof:  If  the  group  is  localizable  z  can  be  decoupled  from  system 

implying  that  n-n  +1  modes  are  unobservable  from  z.  If  V  is  a  basis 

O'  o 

of  the  eigenspace  corresponding  to  these  modes 

G  V  “0  (4.58) 

o  a 


which  implies  that  rows  of  V  corresponding  to  states  in  the  group  are 

o 

equal.  Hence,  it  is  a  coherent  group.  Conversely,  the  rows  of  V 

o 

corresponding  to  a  coherent  group  are  equal  implying  that  n  modes  are 

(T 

unobservable.  When  n  satisfies  (4.57)  the  number  of  observable  modes 

a 

is  n  -1  which  equals  the  dimension  of  z  and  the  group  is  localizable. 

O' 

Note  that  for  V  to  be  full  rank  n  has  to  satisfy 

a  a 


n  <  n  -  n  +  l 
a  ~  a 


(4.59) 


(  The  notion  of  localizability  is  identical  to  aggregability 
with  respect  to  matrix  G0  of  (4.54).  To  avoid  confusion  we  reserve  the 
latter  term  for  "area"  aggregation. 


Hence,  for  the  coherent  group  to  be  localizable  n^  is  required  to  take  its 

maximum  permissible  value.  Note  also  that  the  smaller  the  group  the  larger 

n  and  the  harder  it  is  to  satisfy  the  localizability  conditions.  On  the 
O' 

other  hand  we  do  not  benefit  much  by  localizing  a  large  group.  Finally, 

note  that  the  only  modes  observable  from  the  local  variables  of  a  localiz- 

c 

able  group  are  the  complementary  modes  q  ,  henceforth,  called  local  modes. 

When  a  system  is  divided  into  more  groups  of  states  the  localiz¬ 
ability  conditions  can  be  applied  independently  to  each  group.  As  an 
example,  let  the  modal  matrix  of  a  5-state  system  be 


Xj.  ^2  ^3  *4.  X5 


(4.60) 


m 

» 

®  . 


where  the  stars  can  be  any  numbers  such  that  the  matrix  is  nonsingular. 

Group  a  *  {x^x^x^}  is  a  o^-coherent  group  where  *  {Xj^X^Xj}  and 

satisfies  (4.57).  Hence,  it  is  localizable  and  its  local  modes  are 
c 

°a  “  CX4>X5}.  Likewise  group  g  -  [x4,x5]  is  -  {Xg , X3, X^, X5}  satisfies 

g 

(4.57)  and  it  is  localizable  with  local  modes  er^  *  {X^}.  When,  as  in  the 
above  example  both  the  groups  and  the  local  modes  are  disjoint,  the 
system  is  called  multi- localizable. 
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As  it  is  clear  from  (4.60)  the  multi- localizability  conditions 
are  very  stringent.  Note  that  in  a  multi- localizable  system,  each  set  of 
local  variables  decouples  from  the  rest  of  the  system  and  observes  only 
the  local  modes.  A  less  stringent  requirement  is  that  the  local  variables 
from  the  different  groups  decouple  from  the  system  as  a  single  set. 
Definition  4.9  Let  the  states  of  x=A  be  divided  into  r  disjoint  groups, 
each  consisting  of  two  or  more  states,  and  s  states  not  assigned  to  any 
group.  Then  the  system  is  called  decomposable  if  the  local  variables 


decouple  from  the  system  that  is  if  there  is  a  matrix  A 


4 


such  that 


•R124  422  NONEXPLICIT  SINGULAR  PERTURBATIONS  AND  INTERCONNECTED  2/2 
SVSTENSCU)  ILLINOIS  UNIV  AT  URBANA  DECISION  AND  CONTROL 
LAB  G  H  PEPONIDES  SEP  82  DC-55  N00014-79-C-0424 

f/G  12/1  NL 


UNCLASSIFIED 


:  -  * 


■n  is  the  number  of  states  in  group  v  and  the  s  single  states  are  the 

'v 

last  entries  of  x. 

The  following  lemma  establishes  the  relationship  between 
decomposability  and  coherency. 

Lemma  4.10  A  system  is  decomposable  if  and  only  if  each  of  its  groups  is 
coherent  with  respect  to  the  same  set  of  modes  9 _  and 

cL 

r  +  s  ■  n  (4.64) 

ora 

where  n  the  number  of  modes  in  9  • 
ora  a 

Proof:  If  a  system  is  decomposable,  only  n-s-r  inodes  are  observable  from 

the  same  number  of  local  variables  z.  Let  a  be  the  set  of  r+s  unobservable 

a 

modes  and  V  a  basis  for  the  corresponding  eigenspace.  Then 

GT  V  -  0  (4.65) 

which  implies  that  each  of  the  r  groups  is  a  coherent  group  with  respect 
to  9  .  Conversely  if  all  the  groups  are  coherent  with  respect  to  a  modes 
are  unobservable  from  z.  If  further  the  number  of  the  z  variables  n-s~r 
equals  the  number  of  observable  modes  n-n  ,  that  is  (4.64)  is  satisfied, 
the  system  is  decomposable. 

As  an  Illustration  consider  again  groups  <*  and  3  in  (4.60). 

Both  &  and  0  are  coherent  with  respect  to  oft  ■  *  acf>°£  an<*  t*ie  two 

modes  equal  in  number  the  two  areas.  Hence,  the  system  in  (4.60)  is 


decomposable. 


CHAPTER  5 


SUGGESTIONS  AND  CONCLUSIONS 


5.1  Suggestions  for  Further  Research 

The  ideas  in  this  thesis  can  be  extended  in  several  directions. 

In  Chapters  2  and  3  we  showed  that  equilibrium  and  conservation  properties 
of  an  auxiliary  system  imply  multi-time-scale  behavior.  However,  we  did 
not  investigate  the  relation  between  the  two  properties.  Does  the  existence 
of  one  property  imply  the  other?  And  under  what  conditions?  It  is  clear 
that  in  Linear  Time  Invariant  systems  with  simple  structure  of  7?(Aq) 
(Equations  (2.11)-(2.13))  the  two  properties  are  equivalent.  We  feel  that 
the  existence  of  an  equilibrium  manifold  implies  conservation  properties 
in  a  wide  class  of  nonlinear  systems.  This  issue  and  the  one  of  systematic 
procedures  for  finding  equilibrium  and  dynamic  manifolds  deserve  further 
Investigation.  A  look  at  the  decomposition  in  [63]  and  the  differential 
geometry  techniques  used  therein  should  prove  useful.  Time  scales  in 
discrete-time  systems  is  a  rather  neglected  topic.  Does  the  coordinate- 
free  characterization  carry  over  to  this  class  of  systems?  And  how  are 
time  scales  related  to  the  sampling  period? 

We  have  dealt  mostly  with  time  scales  of  free  systems.  On  the 
other  hand,  high  gain  control  is  known  to  change  the  time-scale  behavior 
(Section  3.3)  of  systems.  In  Section.  2.5  we  have  given  conditions  under 

a 

which  nonexplicit  controlled  LTI  models  can  be  transformed  to  explicit 
controlled  models.  Similar  results  for  nonlinear  systems  would  be  desirable. 


This  issue  is  related  to  extension  of  the  high  gain  results  of  Section  3.3 
to  wider  classes  of  systems .  A  rather  easy  extension  would  be  to  include 
dependence  of  matrix  B  on  x. 

The  decomposition  into  a  slow  core  and  fast  residues  seems 
promising  in  decentralized  and  hierarchical  control  design  along  the  lines 
of  [57].  When  dealing  with  physical  systems  such  as  power  systems, 
prospects  for  Implementation  of  such  designs  should  be  a  consideration. 

In  terms  of  practical  significance  the  time  scale  decomposition 
of  dynamic  networks  in  Sections  4.2,4. 3  seems  to  be  the  most  promising. 
Stability  tests  by  decomposition  methods  have  been  used  in  power  systems 
[64]  but  they  usually  give  conservative  results.  We  feel  that  the  decomposi 
tion  into  slow  core  and  fast  residues  takes  advantage  of  the  structure  of 
the  system  (weakly  connected  areas)  and  it  is  likely  to  provide  practical 
results.  Moreover,  it  can  furnish  information  on  the  type  of  instability, 
that  is  Intermachine  of  interarea  instability.  It  would  also  be  interesting 
to  study  time-scale  separation  and  stability  questions  using  more  complex 
generator  models  and  an  unreduced  network. 


5.2  Conclusions 

Singular  perturbations  have  been  related  to  equilibrium  and 
conservation  properties  of  an  auxiliary  system.  Besides  providing  a 
coordinate-free  characterization  of  singularly  perturbed  systems,  these 
properties  have  been  used  in  definition  of  new  predominantly  slow  and 
predominantly  fast  coordinates.  In  the  new  coordinates  an  extensive  amount 
of  literature  provides  simplified  models,  asymptotic  calculations,  two-stage 
designs  and  stability  tests. 


Results  on  high  gain  feedback  control  have  been  extended  to  a 
class  of  systems  such  larger  than  LT1  systems.  Disturbance  rejection 
behavior  has  been  shown  to  be  similar  to  the  one  In  the  LTI  case. 

lhe  relation  between  weak  connections  and  time  scales  in  a  class 
of  interconnected  systesis  has  been  established.  Weak  connections  combined 
with  equilibrium  and  dynamic  manifolds  of  the  subsystesis  give  rise  to 
multi-time-scale  behavior.  Separation  of  the  time  scales  results  in  a 
slow  core  describing  the  system-wide  dynamics  and  a  set  of  fast  residues 

m 

describing  the  local  dynamics.  In  the  new  representation  recent  stability 
results  can  be  applied  to  give  decentralized  stability  tests. 

When  the  interconnected  system  has  the  added  structure  of  a 
dynamic  network  the  slow  core  and  the  fast  residues  acquire  physical 
significance  and  the  definition  of  slow  and  fast  variables  is  related  to 
physical  laws.  The  linearity  of  these  laws  makes  the  transformation 
separating  the  time  scales  linear.  This  transformation  is  the  area 
aggregation-slow  coherency  one,  developed  for  linearized  models  of  power 
systems. 


APPENDIX  I.  SINGULAR  PERTURBATION  ON  THE  INFINITE  HORIZON 
Under  consideration  is  the  system 

|f  -  f(t,y,z,e) 

(P) 

«  -  g<t,y,r,e) 

with  degenerate  system 

(D)  |f  -  f(t,y,0,0)  ,  y(tQ)  -  yQ 

and  boundary  layer  system 

(BL)  |f  -  g(or»P »z»0) 

where  (or* P)  are  treated  as  parameters.  In  (P)  x,f  are  k-vectors,  g,y 

are  j -vectors  and,  without  loss  of  generality  it  is  assumed  that 

8(t,y,0,0)  ■  0  for  all  t,y. 

Let  |x|  ■  2  Jx^l  be  the  norm  of  x,  let  I  *  [0,»], 
k+1 

SR  ■  {(y»*)  e  E  ;  |y|  +  |*|  <  R}  and  let  sa|y>sR|z  represent  the 

It  1 

restrictions  of  S„  to  E  and  EJ . 

R 

The  following  assumptions  are  made  about  (P),  (D),  (BL). 

(I)  System  (P)  has  a  solution  y  »  y(t),  z  ■  z(t)  that  exists  for 
t_  <  t  <  to. 

o  * 

df4 

(II)  f»8»fv»r  .gj.yg  ,g  €  C  where  f  denotes  the  matrix  r— =■  , 


y<V  '  yo 

2(to>  '  2o 
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(HI)  Function  £  is  continuous  at  z«0,e«0  uniformly  in  (t,y)  €  I  x  SRjy 
and  f(t,y,0,0),  fy(t,y,0,0)  are  bounded  on  I  x  SRjy. 

(IV)  Function  g  is  continuous  at  e»0  uniformly  in  (t,y, z)g  I  x  SR,  and 
g(t,y,z,0)  and  its  derivatives  with  respect  to  t  and  the  components  of  y,z 
are  bounded  on  I  x  SR. 

To  simplify  notation  let  X  be  the  class  of  all  continuous, 
strictly  increasing,  real  valued  functions  d(r),  0  <  r  with  d(0)  «  0; 
and  let  J  be  the  class  of  all  nonnegative,  strictly  decreasing,  continuous, 
real-valued  functions  cr(s),  0  <  s  <  ®  for  which  a(s)-»0  as  s-»®. 

(V)  The  zero  solution  of  D  is  uniform-asymptotically  stable.  That  is, 
if  x  »  } (t,to,yQ)  is  the  solution  of  (D),  3  dg  J  such  that 

•(‘•W  ^  d(ly0l)  c(t"to)  £or  \yo\  *  R» 

(VI)  The  zero  solution  of  (BL)  is  uniform-asymptotically  stable  uniformly 

in  the  parameter  (a,P)€  I  x  SRjy.  i8»  y  “  Y(8»z0»a»P)  is  the 

solution  of  BL,  3  > PS®*  »  such  that 

y(s,zo,Q-,p)  <  e(|zj)  p(s) 

for  all  0  <  s  <  <»,  |yQJ  <  R  and  (g,P)€  I  x  SR  |y. 

Then  the  following  Theorem  is  true  [36]. 

Theorem  [36]  Let  conditions  (I)  through  (VI)  be  satisfied.  Then  for 
sufficiently  small  | y0 |  +  |zq|  and  c  the  solution  of  the  perturbed  system 
(P)  exists  for  tQ  <  t  <  ®,  and  this  Solution  converges  to  the  solution  of 
the  degenerate  system  (D)  as  c-*0+  uniformly  on  all  closed  subsets  of 


t  <  t  <  «. 
o 


APPENDIX  II.  A  STABILITY  THEOREM 


Consider 


(P) 


y  - 

ez  -  g(y,z,e) 


y€  By  c  Rn 

z€  B  c  Rm 
z 


where  B  ,B  denote  closed  spheres  centered  around  y*0,  z*0.  Assume  that 
y  z 

y=0,  z*0  is  the  unique  equilibrium  of  (P)  in  B  *B  and  that  g(y,0,0)  «  0, 

y  z 

for  all  y£  By*  The  reduced  system  of  (P)  is 


(R> 


y  -  f(y,0,0)  A  fr(y) 


and  its  boundary  layer  is 


(BL) 


dz 

dr 


s(y»z(T),o). 


Let  the  following  assumptions  be  satisfied. 

(I)  Reduced  system  (R)  has  a  Lyapunov  function  V:Rn-»R,  such  that  for  all 

y€  By, 

[vy  V(y)]T  fr(y)  <  -  *2(y)  »  >  0 


where  Y(y)  is  a  scalar-value  d  function  of  y  with  y(0)  -  0,  Y(y)  j*  0 
if  y  +  0. 

(II)  Boundary- layer  system  (BL)  has  a  Lyapunov  function  W(y,z):Rnx  Rm-R 
such  that  for  all  y€  B  ,  z6  B 
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%  t 


f’ 


(7Z  W(y,z)]T  g(y,z,0)  <  -  a2  <*2  >  0 

where  C^(z)  is  scalar  valued  and  (f>(0)*0,  (f>(z)^0  if  z#0. 

(Ill)  The  following  inequalities  hold  for  all  y€  B  ,  a€  B 

y  z 

(a)  I?y  W(y,z)]T  f  (y,z)  <  CL  C|)2(z)  +  C2  (J>(z)  Y(y) 

(b)  [V  V(y)]T  [f (y,z)  -  f (y,0)J  <  pL  Y(y)  (J>(z) 

(c)  [vz  W(y,z)]T  [g(y,z,e)  -  g(y,z,0)]  <  eK^  <f>2(z)  +  cK^  Y(y)  (J)  (2) 

Theorem  [20]  If  conditions  (I)- (III)  are  true,  the  origin  x*0,y*>0  is  an 
asymptotically  stable  equilibrium  point  of  (P). 


REFERENCES 


G.  L.  Blankenship,  "Asymptotic  Analysis  in  Mathematical  Physics  and 
Control  Theory,"  in  Special  Issue  on  System  Theory  and  Physics 
(R.  W.  Brockett,  Ed.)»  Ricerche  di  Automatics .  Vol.  X,  pp.  265-315, 
December  1979. 

P.  V.  Kokotovic,  R.  E.  O'Malley,  Jr.,  and  P.  Sannuti,  "Singular 
Perturbations  and  Order  Reduction  in  Control  Theory— An  Overview," 
Automatica.  Vol.  12,  pp.  123-132,  1976. 

J.  H.  Chow  and  P.  V.  Kokotovic,  "A  Two-Stage  Lyapunov-Be liman  Feedback 
Design  of  a  Class  of  Nonlinear  Systems,"  IEEE  Trans,  on  Automatic 
Control.  Vol.  AC-26,  pp.  656-663. 

P.  V.  Kokotovic  and  R.  A.  Yackel,  "Singular  perturbations  of  Linear 
Regulators:  Basic  Theorems,"  IEEE  Trans,  on  Automatic  Control. 

Vol.  AC-17,  pp.  29-37,  1972. 

H.  J.  Kelly,  "Aircraft  Maneuver  Optimization  by  Reduced-Order 
Approximations,"  in  Control  and  Dynamic  Systems,  pp.  131-178, 

Academic  Press,  New  York,  1973. 

F.  G.  Heineken,  H.  M.  Tsuchiya,  and  R.  Aris,  '*)n  the  Mathematical 
Status  of  the  Psuedo-Steady  State  Hypothesis  of  Biochemical  Kinetics," 
Math.  Biosciences.  Vol.  1,  pp.  95-113,  1967. 

D.  S.  Cohen,  Ed.,  Mathematical  Aspects  of  Chemical  and  Biochemical 
Problems  and  Quantum  Chemistry.  SIAM -AMS  Proceedings,  American 
Mathematical  Society,  Providence,  1974. 

R.  E.  O'Malley,  Jr.,  Introduction  to  Singular  Perturbations.  Academic 
Press,  New  York,  1974. 

F.  Hoppensteadt,  "Properties  of  Solutions  of  Ordinary  Differential 
Equations  with  Small  Parameters."  Comm.  Pure  AppI.  Math..  Vol.  XXIV, 
pp.  807-840,  1971. 

V.  F.  Butuzov,  A.  B.  Vasileva,  and  M.  V.  Fedoryuk,  "Asymptotic  Methods 
on  Theory  of  Ordinary  Differential  Equations,"  in  Progress  in  Mathematics 
R.  V.  Gamkrelidze,  Ed.,  Vol.  8,  pp.  1-82,  Plenum  Press,  New  York,  1970. 

A.  B.  Vasileya,  "Singularly  Perturbed  Systems  Containing  Indeterminacy 
in  the  Case  of  Degeneracy."  Soviet  Math.  Dokl..  Vol.  16.  do.  1121-112S 


R.  E.  0  Malley,  Jr.  and  J,  E.  Flaherty,  "Singular  Singular-Perturbation 
Problems,"  Lecture  Notes  in  Mathematics.  No.  594,  Springer-Verlaa 
pp.  422-436.  ** 


13.  S.  L.  Campbell,  Singular  Systems  of  Differential  Equations.  Pitman 
Advanced  Publishing  Program,  1980. 

14.  J.  D.  Cobb,  "Descriptor  Variable  and  Generalized  Singular  Perturbed 
Systems;  a  Geometric  Approach,"  Ph.D.  Dissertation,  University  of 
Illinois,  Urbane,  Illinois,  1980. 

15.  G.  C.  Verghese,  B.  C.  Levy,  and  T.  Kailath,  "A  Generalized  State  Space 
for  Singular  Systems,"  IEEE  Trans,  on  Automatic  Control.  Vol.  AC-26, 
pp.  811-831,  1981. 

16.  N.  Fenichel,  "Geometric  Singular  Perturbation  Theory  for  Ordinary 
Differential  Equations,"  J.  Differential  Equations.  Vol.  31,  pp.  53-98, 
1979. 

17.  Kar-Keung  D.  Young,  Petar  V.  Kokotovic,  and  Vadim  I.  Utkin,  "A  Singular 
Perturbation  Analysis  of  High-Gain  Feedback  Systems,"  IEEE  Trans,  on 
Automatic  Control.  Vol.  AC-22,  No.  6,  December  1977. 

18.  K.  D.  Young,  "Disturbance  Decoupling  by  High  Gain  Feedback,"  to  appear. 

19.  J.  L.  Willems,  "Almost  Invariant  Subspaces;  An  Approach  to  High  Gain 
Feedback  Design,  Part  1:  Almost  Controlled  Invariant  Subspaces," 

IEEE  Trans,  on  Automatic  Control.  Vol.  AC-26,  pp.  235-252,  1981. 

20.  A.  Saber!  and  H.  Khalil,  '^Quadra  tic -Type  Lyapunov  Functions  for 
Singularly  Perturbed  Systems,"  to  appear. 

21.  R.  Fodmore,  "Identification  of  Coherent  Generators  for  Dynamic 
Equivalents,"  IEEE  Trans,  on  Power  Apparatus  and  Systems.  Vol.  PAS-97, 
pp.  1344-1354,  1978. 

22.  F.  Saccomano,  "Dynamic  Modeling  of  Multimachine  Electric  Power 
Systems,"  Svmoosium  on  Mathematical  Methods  for  the  Analysis  of 
Large  Systems.  Praga,  1974. 

• 

23.  K.  N.  Stanton,  "Dynamic  Energy  Balance  Studies  for  Simulation  of 
Power  Frequency  Transients,"  1971  PICA  Proc..  pp.  173-179,  1971. 

24.  B.  Avramovic,  "Time  Scales,  Coherency,  and  Weak  Coupling,"  Ph.D. 
Dissertation,  Report  DC-40,  Coordinated  Science  Laboratory, 

University  of  Illinois,  Urbane,  Illinois,  1980. 

25.  B.  Avramovic,  P.  V.  Kokotovic,  J.  R.  Winkelman,  and  J.  H.  Chow, 

"Area  Decomposition  for  Electromechanical  Models  of  Power  Systems," 
Automatics .  Vol.  16,  pp.  637-648,  1980. 

26.  P.  V.  Kokotovic,  B.  Avramovic,  j.  ...  Ch  ,  and  J.  R.  Winkelman, 
"Coherency  Based  Decomposition  ani-^g rogation,"  to  appear  in 
Automatics.  1982. 


J.  R.  Wlnkelman,  J.  H.  Chow,  B.  C.  Bowler,  B.  Avramovic,  and  P.  V. 
Kokotovic,  "An  Analysis  of  Interarea  Dynamics  of  Multi-Machine 
System,"  IEEE  Trans,  on  Power  Apparatus  and  Systems,  Vol.  PAS-100, 
pp.  754-763,  1981. 

28.  J.  H.  Chow,  B.  Avramovic,  P.  V.  Kokotovic,  and  J.  R.  Winkelman, 
"Singular  Perturbations,  Coherency,  and  Aggregation  of  Dynamic 
Systems,"  Final  Report,  DOE  Contract  DE-AC05-77ET29104,  July  1981. 

29.  G.  Peponides  and  J.  H.  Chow,  "Area  Aggregation  of  Network  Systems," 
Proc.  of  IEEE  Conference  on  Decision  and  Control,  pp.  206-214,  1981. 


C.  A.  Desoer  and  M.  J.  Shensa,  "Networks  with  Very  Small  and  Very 
Large  Parasitics:  Natural  Frequencies  and  Stability,"  Proc .  IEEE , 

Vol.  58,  pp.  1933-1938,  1970. 

S.  S.  Sastry  and  C.  A.  Desoer,  "Jump  Behavior  of  Circuits  and  Systems, 
IEEE  Trans,  on  Circuits  and  Systems,  Vol.  CAS-28,  pp.  1109-1124,  1981. 


P.  A.  Lagerstrom  and  R.G.  Casten,  "Basic  Concepts  Underlying  Singular 
Perturbation  Techniques,"  SIAM  Review.  Vol.  14,  pp.  63-120,  1972. 

L.  A.  Segel,  "Simplification  and  Scaling,"  SIAM  Review.  Vol.  14, 
pp.  547-571,  1972. 

G.  Strang,  Linear  Algebra  and  its  Applications.  Academic  Press, 


New  "  rk,  1976. 

35.  J.  J.  Levin  and  N.  Levinson,  "Singular  Perturbations  of  Nonlinear 
Systems  of  Differential  Equations  and  an  Associated  Boundary  Layer 
Equation,"  J.  Rat.  Mech.  Analyt..  Vol.  3,  pp.  274-280,  1954. 

36.  F.  C.  Hoppensteadt,  "Singular  Perturbations  on  the  Infinite  Interval 
Trans.  Amer.  Math.  Soc..  Vol.  123,  #2,  pp.  521-535,  1966. 

37.  P.  V.  Kokotovic  and  A.  H.  Haddad,  "Controllability  and  Time-Optimal 
Control  of  Systems  with  Slow  and  Fast  Modes,"  Short  Papers,  IEEE 
Trans,  on  Automatic  Control,  February  1975. 

38.  Petar  V.  Kokotovic,  "A  Riccati  Equation  for  Block-Diagonalization 
of  Ill-Conditioned  Systems."  IEEE  Trans,  on  Automatic  Control. 
December  1975. 

39.  Joe  H.  Chow  and  Petar  V.  Kokotovic,  "Two-Time-Scale  Feedback  Design 
of  a  Class  of  Nonlinear  Systems,"  IEEE  Trans,  on  Automatic  Control. 
Vol.  AC-23,  June  1978. 


40.  Hassan  K.  Khalil  and  Petar  V.  Kokotovic,  "Control  of  Linear  Systems 
with  Multiparameter  Singular  Perturbations,"  Antomatica.  Vol.  15, 
pp.  197-207,  Pergamon  Press,  1979. 


99 


41.  K.  W.  Chang,  "Singular  Perturbation  of  a  General  Boundary  Value 
Problem,"  SIAM  J.  Math.  Anal.,  Vol.  3,  pp.  520-526,  1977. 

42.  R.  G.  Phillips  and  P.  V.  Kokotovic,  "A  Singular  Perturbation  Approach 
to  Modeling  and  Control  of  Markov  Chains,"  IEEE  Trans,  on  Automatic 
Control.  Vol.  AC-26,  pp.  1087-1094,  1981. 

43.  P.  J.  Courtois,  Decomposability;  Queuing  and  Computer  System 
Applications .  Academic  Press,  New  York,  1977. 

44.  P.  V.  Kokotovic,  "Subsystems,  Time  Scales,  and  Multimodeling," 
Automatica,  Vol.  17,  pp.  789-795,  November  1981. 

45.  G.  Peponides,  "Two  Time  Scale  Modeling  via  Physical  Considerations," 
Proc.  of  the  15th  Asilomar  Conference  on  Circuits.  Systems  and 
Computers .  1981. 

46.  H.  A.  Simon  and  A.  Ando,  "Aggregation  of  Variables  in  Dynamic  Systems," 
Econometrica.  Vol.  29,  pp.  111-138,  1963. 

47.  A.  N.  Michel  and  R.  K.  Miller,  Qualitative  Analysis  of  Large  Scale 
Dynamical  Systems.  Academic  Press,  New  York,  1977. 

48.  D.  Siljak,  Large  Scale  Dynamic  Systems.  North  Holland,  1978. 

49.  M.  Araki,  "Stability  of  Large  Scale  Nonlinear  Systems -Quadra tic 
Order  Theory  of  Composite  System  Method  Using  M-Matrices," 

IEEE  Trans,  on  Automatic  Control.  Vol.  AC-23,  No.  2,  pp.  129-142, 

1978. 

50.  W.  H.  Fleming,  Functions  of  Several  Variables.  Addison-Wesley, 

Reading,  Massachusetts,  1965. 

51.  L.  H.  Loomis  and  S.  Sternberg,  Advanced  Calculus.  Addison-Wesley, 

Reading,  Massachusetts,  1968. 

52.  U.  Shaked  and  B.  Kouvaritakis ,  "The  Zeros  of  Linear  Optimal  Control 
Systems  and  Their  Role  in  High  Gain  Stability  Design,"  IEEE  Trans, 
on  Automatic  Control.  Vol.  AC-22,  pp.  597-599,  August  1977. 

53.  M.  V.  Meerov,  "Automatic  Control  Systems  Stable  with  Indefinitely 
Large  Values  of  Gain,"  Automatlka  i  Telemekhanika.  Vol.  8,  No.  4, 

1947.  .  . 

* 

» 

54.  U.  Shaked  and  B.  Kouvaritakis,  "Asymptotic  Behavior  of  Root-Locus  of 
Linear  Multivariable  Systems,"  International  Journal  of  Control. 

Vol.  23,  No.  3,  pp.  297-340,  1976^ 

55.  Kar-Keung  David  Young,  "Analysis  and  Synthesis  of  High  Gain  and 
Variable  Structure  Feedback  Systems/'  Ph.D.  Dissertation,  Report  R-800, 
Coordinated  Science  Laboratory,  University  of  Illinois,  Urbana,  Illinois, 
1977. 


56.  N.  R.  Sandell,  Jr.,  P.  Varaiya,  M.  Athans,  and  M.  G.  Safonov, 

"Survey  of  Decentralized  Control  Methods  for  Large  Scale  Systems," 

IEEE  Trans,  on  Automatic  Control.  Vol.  AC-23,  pp.  108-128,  1978. 

57.  Hassan  K.  Khalil  and  Petar  V.  Kokotovic,  "Control  Strategies  for 
Decision  Makers  Using  Different  Models  of  the  Same  System," 

IEEE  Trans,  on  Automatic  Control.  Vol.  AC-23,  No.  2,  April  1978. 

58.  F.  F.  Wu  and  N.  Narsimhamurthi,  "Coherency  Identification  for  Power 
System  Dynamic  Equivalents,"  Memorandum  No.  UCB/ERL  M77/57, 

University  of  California,  Berkeley,  California,  1977. 

59.  S.  Sastry  and  P.  Varaiya,  "Coherency  for  Interconnected  Power  Systems," 
IEEE  Trans,  on  Automatic  Control.  Vol.  AC-26,  pp.  218-226,  1981. 


60.  P.  M.  Anderson  and  A.  A.  Fouad,  Power  System  Control  and  Stability, 

Iowa  State  University  Press,  Ames,  Iowa,  1977. 

61.  J.  H.  Chow,  J.  J.  Allemong,  and  P.  V.  Kokotovic,  "Singular  Perturbation 

Analysis  of  Systems  with  Sustained  High  Frequency  Oscillations," 
Automatica.  Vol.  14,  pp.  271-279,  1978.  j 

62.  J.  J.  Allemong,  "A  Singular  Perturbation  Approach  to  Power  System 
Dynamics,"  Ph.D.  Dissertation,  University  of  Illinois,  Urbana, 

Illinois,  1978. 

63.  A.  Isidori,  A.  J.  Krener,  C.  Gori-Giorgi,  and  S.  Monaco,  "Nonlinear 
Decoupling  via  Feedback:  A  Differential  Geometric  Approach," 

IEEE  Trans,  on  Automatic  Control.  Vol.  AC-26,  pp.  331-345,  1981. 

64.  M.  A.  Pai,  Power  System  Stability.  North  Holland  Publishing  Company, 
Amsterdam,  1981. 

65.  A.  I.  Klimushev  and  N.  N.  Krasovskii,  "Uniform  Asymptotic  Stability  of 
Systems  of  Differential  Equations  With  a  Small  Parameter  in  the 
Derivative  Terms,"  J.  Appl.  Math.  Mech. .  Vol.  25,  pp.  1011-1025, 

1962. 

66.  F.  Hoppensteadt,  "Asymptotic  Stability  in  Singular  Perturbation 
Problems,  II,"  J.  Diff.  Equations.  Vol.  15,  pp.  510-521,  1974. 

67.  R.  E.  O'Malley,  '’On  Two  Methods  of  Solution  for  a  Singularly  Perturbed 
Linear  State  Regulator  Problem,"  SIAM  Review.  Vol.  16,  pp.  16-37, 

1975. 

68.  S.  H.  Javid,  and  P.  V.  Kokotovic,  "A  Decomposition  of  Time  Scales  for 
Iterative  Computation  of  Time-Optimal  Controls,"  Journal  of  Opt. 

Theory  and  Appl..  Vol.  21,  No.  4,  pp.  459-468,  April  1977. 

69.  A.  H.  Haddad,  "Linear  Filtering  of  Singularly  Perturbed  Systems," 

IEEE  Trans,  on  Automatic  Control.  Vol.  AC-31,  pp.  515-519,  1976. 


101 


VITA 

George  Michael  Peponldes  was  born  in  Nicosia,  Cyprus  on  May  16,  1953. 
He  received  Che  Diploma  in  Physics  from  Che  University  of  Athens,  Greece  in 
1976.  While  at  the  University  of  Athens  he  was  supported  by  a  scholarship 
from  the  State  Scholarship  Foundation  of  Greece.  He  has  been  a  graduate 
student  in  the  Department  of  Electrical  Engineering,  University  of  Illinois, 
Urbane  since  fall  1977.  In  January  1979  he  received  the  M.S.  degree  in 
electrical  engineering.  Since  then  he  has  been  a  Ph.D.  candidate  working 
on  multi- time-scale  modeling  and  its  application  to  large  scale  systems  such 
as  power  and  economic  systems .  From  August  1977  to  June  1982  he  has  been 
a  research  assistant  in  the  Decision  and  Control  Laboratory  of  the  Coordinated 
Science  Laboratory. 

Mr.  Peponldes  is  a  member  of  I.E.E.E.  and  the  honorary  society  of 
Phi  Kappa  Phi. 


